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1  Introduction 


Electron  accelerators  capable  of  producing  high-quality  electron  beam  with  energies  of 
10-20  Mega-electron- Volt  (MeV)  are  finding  an  increasing  number  of  applications  in  our 
Society.  Such  applications  include  medicine  (oncology,  and  imaging),  industry  (e.g.  lithog¬ 
raphy),  security,  and  discovery  sciences.  On  another  hand,  over  the  last  two  decade,  needs 
in  elementary  particle  physics  have  driven  the  development  of  new  acceleration  techniques 
capable  of  producing  very  high  accelerating  field  [Giga-electron-Volt  per  meter  (GV/m) 
electric  fields].  These  techniques  could  in  turn  be  transplanted  to  low  and  medium  energy 
accelerator  thereby  resulting  in  tabletop  easy-to-operate  devices  that  could  proliferate  to 
small-scale  laboratories  in  industrial  settings,  hospital  and  Universities.  Beside  enabling  a 
wide  range  of  research  this  table-top  accelerator  could  support  the  accelerator-based  light 
source  capable  of  producing  X  rays. 


The  main  goal  of  the  proposed  research  was  to  explore  the  use  of  wakefield  accelerators 
—  a  concept  developed  in  support  to  high-energy  particle  accelerator  —  as  path  toward  the 
realization  of  compact  accelerator  with  10-20  MeV  output  energy. 
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Figure  1:  Schematic  rendition  of  the  wakefield  driven  X-ray  source  concept.  In  this  version 
of  the  concept  both  passive  bunch  and  wakefield-acelerator  structures  are  dielectric-lined 
waveguides  (DLW’s).  The  red  cylindrical  objects  are  injection  and  ejection  mirrors  for  the 
laser. 

Our  approach  toward  the  development  of  a  compact  beam-driven  accelerator  consists 
of  four  main  components  depicted  in  Fig.  1.  The  production  of  a  bunch  suitable  for  beam- 
driven  acceleration  is  enable  by  the  combination  of  a  conventional  L-band  radiofrequency 
(RF)  photoemission  source  coupled  to  a  “passive  buncher”.  In  this  approach  a  sub-5  MeV 
electron  bunch  emitted  from  the  RF  gun  is  transversely  focused  into  a  high-impedance 
structure  (e.g.  a  dielectric-lined  waveguide,  or  a  corrugated  structure).  The  self-interaction 
of  the  bunch  with  its  short-range  wakefield  imparts  longitudinal-phase-space  correlations 
suitable  for  attaining  high  peak  current  as  the  bunch  drift  in  a  free  space.  The  structure 
geometry  and  electron-beam  parameters  are  optimized  to  form  a  high-peak-current  “drive” 
bunch  followed  by  a  lower  charge  “witness”  bunch  .  The  two-bunch  system  is  then  sent 
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to  a  dielectric-lined  waveguide  (DLW)  with  parameters  optimized  so  that  the  wakefield 
produced  by  the  drive  bunch  accelerates  the  witness  bunch.  In  a  last  stage,  after  disposal 
of  the  drive  bunch,  the  witness  electron  bunch  is  collided  head-on  with  a  ultrafast  laser 
thereby  producing  radiation  with  upshifted  frequency  c o  ~  where  7  is  the  witness 

bunch  Lorentz’  factor  and  cjl  =  2ttc/Xl  with  Xl  being  the  laser  wavelength. 

2  Development  of  efficient  numerical  simulation  of 

WAKEFIELD  ACCELERATION 

2.1  Wakefield  generation 

The  geometry  of  the  problem  analyzed  in  this  paper  is  depicted  in  Fig.  2.  Transient 
effects  resulting  from  the  injection  of  the  electron  bunch  in  the  structure  are  not  included 
(the  structure  is  assumed  to  be  infinitely  long)  and  the  drive  bunch  is  taken  to  be  ultra- 
relativistic  with  its  Lorentz  factor  7  much  greater  than  unity.  The  longitudinal  motion 
of  the  drive  bunch  is  assumed  to  be  at  a  constant  speed  v  and  the  origin  of  the  z- axis  is 
chosen  such  that  z  =  vt  where  t  is  the  time.  The  bunch  is  modelled  has  a  superimposition 
of  line  of  charges  with  finite  extension  in  horizontal  direction  and  with  an  offset  yo  in  the 
vertical  direction. 


p(x,  y,z)  =  \m  cos (kXtmx)S(y  -  y0)S(z  -  vt),  (1) 

m 

The  line-charge  distribution  is  taken  to  be  symmetric  with  respect  to  the  vertical  y- axis 
[p(x,y,z)  —  p(—x,y,z)]  and  vanishes  at  the  dielectric  and  metallic  boundary  surfaces 
[p(— Lx/2,  y,  z)  =  p(Lx/2,  y,  z)  =  p(x,a,z)  =  p(x,—a,z)  =  0].  Under  these  assumptions 
the  charge  distribution  can  be  written  as  a  Fourier  series  in  the  horizontal  direction 

p(x,  y,z)  =  ^2  cos (kX:mx)S(y  -  y0)5(z  -  vt),  (2) 

m 

where  Xm  is  a  constant  and  each  term  is  indexed  by  the  integer  m  —  0, 1,  •  •  •  defined 
such  that  kx  =  (2771  +  1)-^.  In  the  vacuum  region  and  away  outside  the  charge  distribution 
the  wakefield  satisfies  the  wave  equation: 

(v2-^)e  =  0.  (3) 

In  this  paper  we  consider  only  the  propagating  modes  with  harmonic  longitudinal  and 
temporal  dependencies  of  the  form  E  {x,y,z,t)  =  E(x,  y^el(ut-kzz)  wjiere  =  k.z  In 
addition  to  the  propagating-mode  solutions,  Eq.  3  when  subject  to  the  boundary  conditions 
also  admits  evanescent  modes  with  imaginary  eigen-frequencies  uo  .  Since  these  evanescent 
modes  are  present  over  short  distances  behind  the  drive  bunch,  they  do  not  contribute  to  the 
beam  dynamics  of  a  subsequent  witness  bunch.  These  short-range  fields  are  consequently 
ignored  in  the  remaining  of  this  paper. 

In  the  ultra-relativistic  regime  and  considering  the  modes  to  be  synchronous  with  the 
beam  by  imposing  their  phase  velocity  to  satisfy  =  v  ~  c.  The  transverse 
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Figure  2:  Overview  of  the  DLW  structure  geometry  (left)  and  transverse  cross  section 
(right).  A  rectangular-shaped  drive  beam  is  displayed  in  blue. 

wavevector  k±  =  ^k2  —  k2  ~  fc/7  =  0  the  wave  equation  3  can  be  written  only  in  terms 
of  the  partial  derivatives  with  respect  to  the  transverse  coordinates 

Vj_E(a;,  y)  =  0.  (4) 

The  symmetry  of  drive-beam  charge  distribution  and  the  boundary  conditions  at  x  — 
ELX/ 2,  determine  (up  to  constants)  the  analytical  expression  of  the  fields.  In  the  vac¬ 
uum  region.  The  electric  field  components  are  simple  combinations  of  trigonometric  and 
hyperbolic  functions: 


Ex 

OC 

sin  (kX}mx) 

cosh  (kyy), 

Ey 

oc 

cos  (kx,mx) 

Isinh  (kyy), 

(5) 

Ez 

oc 

cos  (kx,mx) 

1  cosh  (kyy). 

The  axial  field  Ez  associated  with  this  set  of  solution  is  symmetric  with  respect  to  the  hor¬ 
izontal  axis  [i.e.  Ez(x,—y,z)  —  Ez{x,y,z)\  we  henceforth  refer  this  set  to  as  “monopole” 
modes.  Similarly,  a  set  of  field  with  an  antisymmetric  axial  field  [. Ez(x ,  — y,  z)  =  —Ez(x,  y,  z)] 
is  obtained  by  substitution  of  sinh (kyy)  with  cosh (kyy)  and  vice  versa.  This  latter  set  is 
termed  as  “dipole”  modes  in  the  remaining  of  this  paper. 

Inside  the  dielectric,  the  transverse  wave  number  cannot  be  neglected 

cj2 

k±  =  kl,m  +  kl~  !)  >  0.  (6) 

where  er  is  the  relative  electric  permittivity  of  the  dielectric  medium.  The  boundary 
conditions  at  y  —  Eb  and  at  x  =  ELx/2  determine  the  trigonometric  form  of  the  field 
expressions.  Inside  the  dielectric  region,  there  is  no  distinction  between  monopole  and 
dipole  modes  and  the  electric  field  components  are  given  by 

Ex  oc  sin(kXjmx)  cos[ky(b  -  y)}, 

Ey  oc  cos(kXjmx)sin[ky(b  -  y)}, 

Ez  oc  cos(kXjmx)cos[ky(b-y)]. 

The  corresponding  expressions  for  the  magnetic  field  in  the  vacuum  and  dielectric  regions 
can  be  easily  obtained  following  a  similar  prescription. 


(7) 
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Inspection  of  the  z-component  of  the  fields  indicates  that  the  normal  modes  cannot 
be  categorized  in  the  usual  transverse  electric  or  transverse  magnetic  sets.  The  reason  is 
that  the  separation  surface  between  the  dielectric  and  vacuum  regions  is  in  the  x  —  z  plane 
unlike  the  case  of  the  uniformly-filled  waveguides  where  this  surface  is  in  the  transverse 
x  —  y  plane. 

Therefore,  it  is  natural  to  categorize  the  modes  depending  on  the  orientation  of  the  fields 
with  respect  to  the  dielectric  surface.  Following  the  definition  introduced  in  Ref.  we  clas¬ 
sify  the  mode  as  Longitudinal  Section  Magnetic  (LSM)  and  Longitudinal  Section  Electric 
(LSE)  modes  corresponding  respectively  to  the  case  when  the  magnetic  and  electric  field 
component  perpendicular  to  the  dielectric  surface  is  vanishes.  In  the  configuration  shown 
in  Fig.  2,  LSE  and  LSM  modes  correspond  respectively  to  Ey  —  0  and  Hy  =  0  at  the 
vacuum-dielectric  interface. 

2.2  Dispersion  equation  &  electromagnetic  Wakefields 
2.2.1  Longitudinal  Section  Magnetic  (LSM)  modes 

To  obtain  the  normal  mode  frequencies  it  is  convenient  to  use  the  hertzian  potential  vector 
method  .  In  the  source-free  region,  magnetic  and  electric  fields  can  be  expressed  in  terms 
of  the  vector  potential  II: 


H  =  icceV  x  II,  (8) 

e  =  fc2n  +  v(v-n),  (9) 

where  II  is  the  solution  of  the  homogeneous  Helmholtz  equation  in  the  vacuum  region. 
To  determine  the  frequencies  and  the  fields  for  the  modes  when  Hy  —  0  (LSM)  it  is 
convenient  to  choose  II  =  IIe  =  ripe(x:y)el^ujt~kzZ")y  where  the  the  scalar  function  ^e{x^y) 
must  be  a  solution  of  (Vj_  +  =  0.  From  Eqs.  8  and  9  all  fields  can  be  expressed  in 

terms  of  the  unknown  function  ^e(x,y): 


rp  _  d2ipe 

x  dxdy 

Ey  =  fcVe  +  ^T 

Ez  =  - ikj ^ 


Hx  =  —ek2zcipe 
dip, 


Hy=  0 


(10) 


Hz  —  iekzc 

ox 


The  expression  for  ipe  that  satisfies  Eqns.  10,  6,  and  8  is  given  by: 


Ipe  = 


Acos(kx,mx)  sinh (kXjTny), 

B  cos (kX;mx)  cos  [ky(b  -  y)\ , 


0  <  y  <  a 
a  <  y  <  b 


(ii) 


where  A  and  B  are  constants.  The  boundary  conditions  completely  determine  the 
eigen-frequencies  of  the  LSM  modes.  Since  Hx  (oc  e^e)  and  Ez  (oc  ^^)  are  continuous  at 
y  =  a  the  constants  A  and  B  can  be  eliminated  giving  the  dispersion  equation 
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coth (kx^ma)  cot  [ky(b  —  a)]  — 


trkx.m 


(12) 


Therefore  each  discrete  value  of  kx,m  there  is  an  infinite  set  of  discrete  ky,n  values  where 
n  is  an  integer.  So,  the  eigenfrequencies  are  indexed  by  the  integer  couple  (m,  n)  and  verify 


hz  4.  hz  = 

.x.m  1  '^y^n 


(jO, 


m,n 


(6r  !)• 


(13) 


It  worthwhile  noting  that  the  boundary  conditions  do  not  completely  determine  either 
the  function  ri\)e  or  the  fields.  The  reason  being  that,  up  to  this  point,  the  field  source 
terms  (p,j)  were  not  taken  into  account  although  their  symmetries  were  invoked.  Still, 
it  is  straightforward  to  show  that  all  fields  associated  to  on  set  of  mode  (and  also  jpe) 
depend  have  a  common  normalization  constant,  the  amplitude  i?o;m,n5  which  remains  to 
be  determined.  The  expressions  of  the  fields  are  given  by 


E 


x,m,n 


E, 


y,m,n 


Ez 


Hr 


H 


0  <  y  <  a 
a  <  y  <  b 


f  _iEo-,mg.kx,m  gin^^x)  cosh (kx>my) 

MKmx)  sin  [ky,n{b  -  y)\ 

EeC  (kx,m  +  kl)  cos  (kx,mx)  sinh  (kx>my)  0<y<  a 

jSeTs  ^  +  klm)  cos  {kx,mx)  cos  [ky,n(b  -  y)]  a<y<b 

(  E0-m,n  cos (kx>mx)  cosh(kx,my) 

1  ;m,n sin [ky,n(b-a)]  cos (kx,mx)  sin  \ky,n{b  —  y)\  a  <  y 

[  cos(kx,mx)  sinh (kX:Tny)  0  <  y  <  a 

l  - &Sl  cos (kx,mx)  cos  [ky,n{b  -y)\  a  <  y  <  b 


0  <  y  <  a 
f  <  b 


y,m,n 


Hz 


Eq 

;m,n  ec  sm(kXtTnx)  sinh(kx,my) 

Eo;m,nkx,m£C  COsh  {kx,m^) 


0  <  y  <  a 

Fwb-a)]  sin(kx>mx)  cos  [kytTl(b  -  y)\  a  <  y  <  b 


(14) 


2.2.2  Longitudinal  Section  Electric  (LSE)  modes 

The  case  of  the  LSE  modes  ( Ey  =  0)  can  be  treated  in  a  similar  fashin  than  for  the  LSM 
modes.  The  hertzian  vector  electric  potential  is  replaced  by  a  vector  magnetic  potential 
Uy  which  is  related  to  the  fields: 


E  =  — icj/iV  x  IIh 

h  =  k2nh  +  v(v  •  nh) 


(15) 

(16) 
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As  in  the  previous  case,  it  is  convenient  to  factor  out  the  t  and  ^-dependencies  of  the 
Ilh  and  to  define  a  scalar  function  ^(x,  y):  Ilh  =  ^h(x^  y)eL(ujt~kzZ^y.  It  is  straightforward 
to  derive  the  dispersion  equation  for  the  LSE  modes 

k 

coth (kxa)  cot  [ky(b  —  a)]  =  —  (17) 

ky 

and  the  expressions  for  the  electromagnetic-field  components: 


Ex 


lE°C,mkz  s'm(kx,mx)  cosh (kX)m,y) 

Si»(  W)  si“  hn(t>  -  #)] 


0  <  y  <  a 
a  <  y  <  b 


E, 


y,m,n 


o 


Ez 


H, 


x,m,n 


Hi 


y,m,n 


Hz 


.  '  '  0111  [i^y  }nu  u,)\ 

cos(kx,mx)  smh(kx>my) 

coS(fe,raI)  cos  [fc„,„(6  -  y)] 

0  <  y  <  a 


kzflC  - v~u;,y /*-/  ~ 

cosh(fea;;ma) 

kzflC  sin[fcy,n  (6-0)] 

sin<A,™x)  c°Mkx,my) 

(J4  ,m  +fcz)  COSh^/Ca^  777,(2)  •  /7  x 

kx,mkzfj>c  sin [ky,n(b-a)\  b  v  x’m' 


a  <  y 


nXJ 


cos  [ky>n(b  -  y)\ 


kx,mkz  fAC 

sin (kx,mx)  smh.(kX}Tny) 

.  -  sSK^-a)]  cos  [  (b  -  y)\ 


a  <  y  <  b 


0  <  y  <  a 
a  <  y  <  b 


(18) 


2.2.3  Discussion 

For  higher-order  modes  the  solutions  of  both  LSM  and  LSE  dispersion  equations  tend  to 
be  equally  spaced.  For  large  values  of  x  =  ky(b  —  a),  the  solutions  of  the  transcendental 
equations  12  and  17  approach  pn  and  (p  +  ^)  n  respectively  where  the  p^>  1  is  an  integer. 

The  dispersion  equations  for  the  LSM  and  LSE  dipole  modes  is  obtained  by  substituting 
tanh (kx^ma)  for  coth (kx^ma)  in  Equations.  12  and  17  respectively. 

An  important  property  of  the  wakefields  described  by  Eqns.  14  and  18  is  that  the  phase 
difference  between  the  longitudinal  and  transverse  components  is  always  ±|  no  matter 
what  the  mode  type  is.  This  means  that  whenever  the  longitudinal  force,  Fz  =  qEZ:  reaches 
a  minimum  (drive  bunch)  or  a  maximum  (test  bunch),  the  vertical  force,  Fy  —  q  (Ey  +  vBx ) 
is  zero.  The  important  consequence  is  that  the  beam  breakup  cannot  happen  for  these  kind 
of  structures  provided  the  electron  bunch  length  is  much  shorter  than  the  wavelength  of 
the  fundamental  mode. 

2.3  Wakefield  amplitudes 

Although  we  have  obtained  the  general  expression  for  the  electromagnetic  field,  the  con¬ 
stant  Eo;m,ri  still  remained  to  be  determined.  An  often  used  method  to  find  the  constant 
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Eo;m,n  in  Eqns.  14  and  18,  is  to  evaluate  the  Green  function  of  the  wave  equation  with 
sources  included.  The  latter  is  usually  a  cumbersome  procedure  and  we  instead  chose 
to  determine  based  on  energy  balance  considerations.  The  electromagnetic  energy 

stored  in  the  DLW  equals  the  mechanical  work  performed  on  the  drive  bunch.  The  key 
point  here  is  to  invoke  the  fundamental  wakefield  theorem  .  For  a  linear  wakefield,  the 
decelerating  field  acting  on  a  drive  point-charge  is  half  of  the  wakefield  amplitude,  i.e. 
Ed  =  E°’™’n .  Suppose  the  drive  charge  moves  over  an  infinitely  short  distance  Sz.  The 
work  performed  by  the  wakefield  on  the  drive  charge  should  equate  the  energy  stored  in 
the  field 


E  j  - —dxdySz 

m,n  ~~ 

=  ^  J  {eE2  +  /i0 H2)  dxdy5z ,  (19) 

where  the  integration  extends  over  the  transverse  plane.  It  is  straightforward  to  evaluate 
the  full  expressions  of  the  wakefield  amplitudes  in  terms  of  the  wave  numbers  kx,m  and 
ky:U  from  Eqns.  19,  14  and  18: 


t?LSM 
^0  ;m,n 


l 


Am  cosh(fca,jm?/o) 


2eo  sinh(2fca;,ma) 


+ 


€r  cosh2(fea;;yna) 
sin  2[ky,n(p-a)\ 


sin[2  ky,n{b—a)] 
^ky^n 


U2 

^y,  n 


)} 


J^LSE 
^0  ;ra,n 


cosh(kx,my0) 


2eo  sinh(2fe  g ;,ma) 
2  kx,m 


+ 


coslV 
sin2  [h 


{kx,ma)  f  b — a  (  I  k; y , n  \  sin[2 ky^n(b  a)] 

j,n(b— a)\  \  2  yr  /  ^ky,n 


(20) 


The  field  amplitudes  for  the  dipole  modes  can  be  obtained  from  the  previous  equations 
by  substituting  sinh (kX:Tnu)  for  cosh^^Ti)  where  u  takes  on  the  values  of  a  and  y®. 
As  expected,  the  wakefield  amplitude  scales  linearly  with  the  drive  bunch  charge  and 
inverse  proportionally  with  the  transverse  size  of  the  structure.  To  this  point  this  model  is 
fully  three-dimensional  and  the  only  limitation  stems  from  the  assumed  symmetry  of  the 
transverse  drive  charge  distribution  with  respect  to  the  y  axis. 

2.4  Comparison  with  three-dimensional  FDTD  simulations 

To  evaluate  the  wakefield  associated  to  a  drive  bunch,  an  integration  over  the  full  three- 
dimensional  continuous  charge  distribution  must  be  performed.  In  practice  the  integration 
is  replaced  by  numerical  summations  of  discrete  charge  distributions  similar  to  those  de¬ 
scribed  by  Eq.  2.  This  process  is  similar  to  the  charge  discretization  procedure  used  in 
standard  PIC  algorithms  .  An  important  feature  of  this  model  is  that  the  integration  over 
x-direction  is  already  performed  through  the  Fourier  expansion  of  the  charge  distribution. 
For  the  charge  distributions  with  the  considered  in  this  paper,  only  a  few  Fourier  terms 
(<  10)  are  sufficient  to  obtain  an  accurate  representation  of  the  distribution  along  the  x 
direction.  This  is  significantly  less  than  the  number  of  grid  points  in  x-direction  needed 
by  most  PIC  codes  to  evaluate  the  3D-collective  effects  and  external  fields. 
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Figure  3:  Longitudinal  (a)  and  transverse  sections  (b)  of  the  DLW  model  implemented 
in  VORPAL.  The  rectangular  box  delimited  by  green  blocks  represents  the  3D  computa¬ 
tional  domain  used  in  the  simulations.  Grey,  green  and  purple  blocks  respectively  stand 
for  perfectly-matched  layer  (PML),  perfectly-conducting  boundary  (PCB),  and  dielectric 
material  (with  associated  relative  dielectric  permittivity  er).  The  cyan  rectangles  represent 
the  electron  bunch  distribution. 

The  integration  over  the  vertical  direction  is  straightforward  and  in  most  cases,  when 
the  drive  charge  distribution  is  also  symmetric  with  respect  to  the  horizontal  axis,  the 
contribution  of  the  dipole  modes  cancels  out. 

Since  the  phase  velocity  of  the  wakefield  is  the  same  as  the  velocity  of  the  drive  beam, 
causality  principle  requires  that  the  wakefield  vanishes  ahead  of  the  drive  charge.  Therefore, 
the  integration  over  the  longitudinal  direction  extends  only  from  the  observation  point  to 
the  actual  drive  charge  position.  The  wakefield  assumes  the  form 

/oo 

p(z')Wm,n{z  -  z')dz' ,  (21) 

m=U,±,"-  74=u,l,-”  •' 

where  W  stands  for  any  of  the  electromagnetic  field  components  and  Wm^n  are  the  corre¬ 
sponding  field  component  associated  to  the  LSMm:n  and  LSEm5n  modes  given  respectively 
by  Eq.  14  and  Eq.  18.  The  summation  is  performed  over  all  modes  excited  modes. 

The  results  obtained  from  Eq.  21  are  benchmarked  against  simulations  performed  with 
VORPAL  a  conformal  FDTD  (CFDTD)  PIC  electromagnetic  solver  .  vorpal  is  a  parallel, 
object-oriented  framework  for  three  dimensional  relativistic  electrostatic  and  electromag¬ 
netic  plasma  simulation.  The  DLW  model  implemented  in  VORPAL  is  fully  three  dimen¬ 
sional;  see  Fig.  3.  The  model  consists  of  the  rectangular  DLW  surrounded  by  perfectly- 
conducting  boundaries  (PCBs),  The  lower  and  upper  z  planes  are  terminated  by  perfectly 
matched  layers  (PMLs)  that  significantly  suppresses  artificial  reflections  of  incident  radia¬ 
tion  .  A  particle  source  located  on  the  surface  of  the  lower  z  plane  produced  macroparticles 
uniformly  distributed  in  the  transverse  plane  (x,  y)  and  following  a  Gaussian  longitudinal 
distribution.  The  longitudinal  Gaussian  distribution  is  truncated  at  ±3 az  where  <jz  is 
its  root-mean-square  (rms)  length.  Finally,  a  “particle  sink”  at  the  upper  z  plane  allows 
macroparticles  to  exit  the  computational  domain  without  being  scattered  or  creating  other 
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Figure  4:  Eigenfrequencies  [fm^n  =  ujm:n/ (2n)\  and  associated  amplitudes  for  the  mode 
induced  in  the  DLW  with  parameter  listed  in  Tab.  1.  The  (m,  n)  LSM  and  LSE  modes  are 
respectively  shown  as  red  and  blue  bars  . 


source  of  radiation. 

In  order  to  precisely  benchmark  our  model,  a  DLW  which  support  LSM  and  LSE  modes 
is  chosen.  The  parameters  of  the  structure  and  driving  bunch  are  gathered  in  Tab.  1  and 
the  frequency  and  amplitude  associated  to  the  first  few  modes  appear  in  Fig.  4.  The  drive- 
bunch  energy  is  set  to  £  =  1.0  GeV  consistently  with  the  ultra-relativistic  approximation 
used  in  the  analytical  model. 

The  longitudinal  component  of  the  electric  field  simulated  with  VORPAL  is  shown  in 
Fig.  5  as  a  two-dimensional  projection  in  the  y  —  z  plane. 

The  comparison  between  theoretical  calculation  and  VORPAL  simulations  are  shown  in 
Figs.  6  and  7.  The  fields  in  Fig.  6  are  evaluated  as  a  function  of  the  axial  coordinate 
z  at  a  given  transverse  location  (x  —  3.0  mm,  y  —  2.0  mm)  which  corresponds  to  the 
upper-left  corner  of  the  charge  distribution.  Figure  7  displays  the  electromagnetic  field 
evaluated  as  a  function  of  the  horizontal  [Fig.  7  (left)]  and  vertical  [Fig.  7  (right)]  transverse 
coordinate  at  a  given  axial  location  z  =  18.6  mm  corresponding  to  the  minimum  Ez  shown 
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Figure  5:  Snapshot  of  the  axial  electric  field  Ez(x  =  0,7/,  z)  in  the  mid-plane  of  a  slab 
DLW  for  the  bunch  and  structure  parameters  shown  in  Tab.  1.  The  field  was  obtained 
from  VORPAL  simulations. 

in  Fig.  6  (i.e.  maximum  accelerating  field),  fixed.  All  plots  show  a  decent  agreement 
(relative  discrepancy  <  25  %)  between  our  model  and  VORPAL  simulations.  The  notable 
disagreement  observed  in  Fig.  6  for  the  transverse  fields  in  the  vicinity  of  the  driving  charge 
(i.e.  z  27.5  mm)  is  rooted  in  the  absence  of  velocity  fields  in  our  model  (only  radiation 
field  contributes  to  the  wakefield)  while  the  VORPAL  simulations  include  both  velocity 
and  radiation  fields.  In  fact  given  the  bunch  charge  and  duration,  the  amplitude  of  the 
velocity  field  can  be  evaluated  by  convolving  the  charge  distribution  with  the  electric  field 
generated  by  an  ultra-relativistic  particle  E(r,  Q  —  2g(r^/rjJ5(£)  where  q  is  the  particle’s 
charge  and  =  (x,y)  and  (  are  respectively  the  transverse  and  longitudinal  coordinates 
of  the  observation  point  referenced  to  the  particle’s  location.  Such  a  convolution  with  our 
transversely-uniform  longitudinally-Gaussian  charge  distribution  yields  the  electric  field  at 
(x,7/,£)  =  (wx,wy,  0)  to  be  resulting  in  a  value  of  ^  0.8  MV/m  for  the  bunch  parameter 
listed  in  Tab.  1.  The  latter  value  is  in  agreement  with  the  observed  difference  between  the 
VORPAL  and  theoretical  models;  see  Ex  and  Ey  components  in  Fig.  7. 

2.5  Two-dimensional  limit 

An  interesting  limiting  case  occurs  when  Lx  Ly.  In  this  “two-dimensional  limit”,  the 
dependence  of  the  fields  on  the  horizontal  coordinate  x  is  weak  and  completely  vanishes 
when  Lx  — »  oo  (so  that  kx  ~  0  and  mn  —  0).  In  such  a  case  the  field  amplitudes  associated 
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Table  1:  Parameter  of  the  DLW  structure  and  drive  bunch  used  for  benchmarking  our 
theoretical  model  with  VORPAL  simulations. 


parameter 

symbol 

value 

unit 

vacuum  gap 

a 

2.5 

mm 

height 

b 

5.0 

mm 

width 

Lx 

10.0 

mm 

relative  permittivity 

6r 

4.0 

— 

rms  bunch  length 

°z 

1.0 

mm 

bunch  charge 

Q 

1.0 

nC 

bunch  full  width 

Wx 

6.0 

mm 

bunch  full  height 

Wy 

4.0 

mm 

bunch  energy 

£ 

1 

GeV 

to  the  LSM  and  LSE  modes  are  respectively 


T^LSM 

^0;0,n 


4A 


a  + 


er(b—a)  ’ 
sin2  [ky{b—a)\ 


and  E 


LSE 

0;0,n 


~  0, 


(22) 


where  A  is  charge  per  unit  length  in  the  horizontal  direction.  The  LSE  modes  are 
suppressed  and  the  latter  equation  is  in  agreement  with  the  results  of  Ref.  . 

This  limit  case  can  be  practically  reached  by  using  structures  with  large  aspect  ratios 
( Lx  Ly)  driven  by  flat  electron  beams  (tailored  such  that  ax  ay ).  Flat  beams  can 
be  produced  in  photoinjectors  by  using  a  round-to-flat  beam  transformation  .  In  such 
a  scheme,  a  beam  with  large  angular-momentum  is  produced  in  a  photoinjector  .  Upon 
removal  of  the  angular  momentum  by  applying  a  torque  with  a  set  of  skew  quadrupole, 
the  beam  has  its  transverse  emittance  repartitioned  with  a  tunable  transverse  emittances 
ratio  .  Flat  beams  with  transverse  sizes  of  ~  100  fim  and  aspect  ratio  of  ^  20  have  been 
produced,  even  at  a  relatively-low  energy  of  15  MeV  . 

Based  on  the  previous  experience  in  producing  flat  beams  and  preliminary  simulation 
of  the  Advanced  Superconducting  Test  Accelerator  (ASTA)  currently  in  construction  at 
Fermilab.  It  is  reasonable  to  consider  a  250-MeV  3-nC  flat  beam  generated  from  a  pho¬ 
toinjector  to  have  a  the  following  parameters  ax  =  XX,  ay  =  40  /im,  az  =  50  /i m. 
Considering  a  structure  with  a  —  100  /im,  b  —  300  /im  and  e  =  4.5  would  yield  a  maximum 
axial  wakefield  amplitude  of  ^  300  MV/m;  see  Fig.  8.  The  fundamental  LSM  frequency 
is  /i,i  =  193  GHz  and  the  corresponding  LSE  mode  amplitude  is  ^  136  times  lower.  In 
Fig.  8  the  wakefield  is  computed  with  the  asymptotic  limit  provided  in  Eq.  22  and  yet  is 
in  excellent  agreement  with  the  FDTD  simulations. 

2.6  IMPLEMENTATION  IN  A  PARTICLE-IN-CELL  BEAM  DYNAMICS  PRO¬ 
GRAM 


Although  FDTD  simulations  provide  important  insight  that  can  aid  the  design  and  opti¬ 
mization  of  the  DLW  geometry,  their  use  to  optimize  a  whole  linear  accelerator  would  be 
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Figure  6:  Comparison  of  the  electromagnetic  field  components  calculated  with  our  model 
(blue  dashed  line)  and  simulated  with  VORPAL  (red  solid  lines).  The  fields  are  computed 
on  a  line  parallel  to  the  2:  axis  with  transverse  offset  x  =  3.0  mm  and  y  =  2.0  mm.  The 
center  of  the  drive  beam  is  at  2:  =  27.5  mm.  The  DLW  and  bunch  parameters  are  the  one 
displayed  in  Tab.  1. 


time  and  CPU  prohibitive.  Therefore  it  is  worthwhile  to  include  a  semi-analytical  version 
of  the  model  developed  in  this  paper  in  a  well-established  beam  dynamics  program  IMPACT- 
T  .  Our  main  motivation  toward  this  choice  stems  from  the  availability  of  a  wide  range  of 
beamline  elements  models.  In  addition,  IMPACT-T  takes  into  account  space-charge  forces 
using  a  three-dimensional  electrostatic  solver.  The  algorithm  consists  in  solving  Poisson’s 
equation  in  the  bunch’s  rest  frame  and  Lorentz-boosting  the  computed  electrostatic  fields 
in  the  laboratory  frame. 

In  a  typical  particle-tracking  PIC  program,  like  IMPACT-T,  an  electron  bunch  is  de¬ 
scribed  by  a  set  of  “macroparticles”  with  arranged  to  mimic  the  bunch  phase  space  dis¬ 
tribution.  Each  macroparticle  represents  a  large  number  of  electrons  (typically  103  in  our 
simulations).  To  evaluate  the  electrostatic  fields  in  the  bunch’s  rest  frame  the  macroparti¬ 
cles  are  deposited  on  the  cells  of  a  three-dimensional  grid. 

In  order  to  implement  the  model  described  in  this  paper  it  is  convenient  to  replace  the 
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Figure  7:  Comparison  of  the  electromagnetic  field  components  calculated  with  our  model 
(blue  dashed  line)  and  simulated  with  VORPAL  (red  solid  lines).  The  fields  are  computed  a 
a  given  longitudinal  position  z  =  18.6  mm.  The  left-column  plots  display  the  dependence 
on  the  horizontal  coordinate  x  at  y  —  0.043  mm  while  the  right-column  plots  show  the 
dependence  on  the  vertical  coordinate  y  at  x  —  0.043  mm.  The  DLW  and  bunch  parameters 
are  the  one  displayed  in  Tab.  1. 


Fourier  expansion  coefficients  Am  from  Eq.  2  with  discrete  charges  at  the  grid  nodes.  Con¬ 
sidering  a  charge  distribution  composed  two  symmetrically  located  point-charges,  p{pc)  — 
q  [6(x  —  xq)  +  5(x  +  xq)]  then 


A 


m 


2q\ 

XO  COsiJc^rnX  q) 


L 


X 


(23) 


In  general,  the  Am’s  are  sums  over  the  nodes  with  x  >  0.  In  practice  the  sums  extend 
over  all  nodes  provided  that  the  charges  located  at  symmetric  nodes  are  about  the  same 
Qi(xo)  ~  q2(-x0). 

The  general  expressions  of  the  wakefields  at  a  given  position  have  the  following  form 
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Figure  8:  Longitudinal  electric  field  for  a  DLW  with  Lx  Ly  driven  by  a  3-nC  bunch.  The 
bunch  has  a  rectangular  transverse  shape  with  full  width  wx  =  3.0  mm  and  wy  —  150  /im, 
and  a  Gaussian  longitudinal  distribution  with  rms  length  crz  =  50  /im.  The  theoretical 
model  based  on  Eq.  22  (dashed  blue  line)  is  compared  with  VORPAL  simulations  (red  line). 


Ei  =  fi(x,  y,  n)  (  cos  [kz(z0  -  z)]Ti(z0;n )  j  ,  (24) 

n  \zo>z  J 

where  i  =  x,y,z:  n  is  the  mode  index,  fi  s  are  geometry  dependent  functions  of  the  type 
described  in  Eq.  6  and 


Ti(z0;n )  =  E  Q\x0,yo,zo  cos(kx,mx0)cosh(kx,my0).  (25) 

xo,yo 

In  the  latter  equation  the  substitution  cosh  sinh  should  be  made  when  considering 
dipole  modes. 

The  simulated  electromagnetic  field  components  are  compared  with  VORPAL  simula¬ 
tions  in  Fig.  9  for  the  same  case  as  presented  in  Fig.  6.  In  the  Fig.  9  the  IMPACT-T 
simulation  are  performed  with  and  without  activating  the  space-charge  algorithm.  When 
accounting  for  space-charge  forces,  IMPACT-T  is  in  very  good  agreement  with  VORPAL. 

The  convolution  over  the  ^-summation  in  Eq.  24  can  be  be  performed  with  a  numerical 
Fast  Fourier  Transformation  (FFT).  Therefore  the  total  number  of  operations  needed  to 
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Figure  9:  Comparison  of  the  electromagnetic  field  components  obtained  from  IMPACT  T 
[with  (green  dashed  line)  and  without  (blue  dashed  line)  accounting  for  space  charge  forces] 
with  the  fields  simulated  with  VORPAL  (red  solid  lines).  The  fields  are  computed  on  a  line 
parallel  to  the  z  axis  with  transverse  offset  x  =  3.0  mm  and  y  —  2.0  mm.  The  center  of 
the  drive  beam  is  at  z  =  27.5  mm.  The  DLW  and  bunch  parameters  are  the  one  displayed 
in  Tab.  1. 

evaluate  the  wakefields  is  oc  Nmo(iesN^NyNz  log  Nz.  The  total  number  of  modes  is  twice 
the  product  between  the  modes  allowed  for  each  of  the  transverse  wave  numbers:  Nmodes  = 
2NkxNky.  The  factor  of  2  comes  from  the  inclusion  of  the  dipole  modes  along  with  the 
vertically  symmetric  monopole  modes.  For  the  data  generated  in  Fig.  9,  the  IMPACT-T 
simulations  are  more  than  two  orders  of  magnitude  faster  than  the  VORPAL  ones. 
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3  Electron-bunch  shaping  for  enhanced  transformer 
ratio 


3.1  Tradeoffs  between  transformer  ratio  and  peak  acceler¬ 
ating  FIELD 


The  maximization  of  1Z  has  been  well  studied  and  can  be  achieved  by  making  the  decel¬ 
erating  field  constant  over  the  drive  bunch.  In  contrast,  the  maximization  of  E+  can  be 
achieved  by  maximizing  the  peak  current.  This  however,  produces  a  very  large  E_  which 
effectively  reduces  1Z.  We  consequently  investigated  the  relationship  between  E+  and  1Z 
for  different  current  shapes.  The  six  shapes  considered  are  as  follows: 


Gaussian  distribution:  g{x)  =  ^  ^ ^2(j2\  where  g  and  a  are  the  mean  and 

root-mean-square  (rms)  length, 


•  Linearly-ramped:  g[z )  =  mz  +  6,  with  m  and  b  as  parameters, 

5 

•  Fourier  distribution:  g(z )  =  bn  sin(n2:)  where  bn’s  are  free  parameters, 

71=1 

•  Double  triangle  distribution  as  a  combination  of  two  linearly-ramped  distribution 
parametrized  by  ln  and  hni  the  horizontal  and  vertical  position  of  the  n-th  vertice 
respectively. 


•  Exponential  distribution:  g[z)  —  e  az  for  z  G  [0,  L\  with  a  and  L  as  parameters, 

•  Skewed  gaussian  distribution:  g(z)  =  where  £  is  a  shift,  c o 

_El 

the  characteristic  length,  and  a  the  skew,  <fi(x)  =  ^=e  2  ?  and  <fr(x)  =  1 


1  + erf  ( % 


■‘Realistic’1  current  achievable  with  longitudinal-phase-space  manipulation  g(z)  = 


(o/+A1/2(zj))2 


2 - 1  x  @[A(z)l,  with 


fdziy(z)ex where  IJ(z)  =  ^ly^expt - ^ 

A (z)  =  a2  +  4:bfZ  and  ©()  is  the  Heaviside  function.  The  final  current  shape  is 
therefore  controlled  via  the  parameters  aj  and  bf.  Here  we  take  au  —  0.05. 


Typical  shapes  associated  to  these  functions  g(z)  are  displayed  in  Fig.  10  (left).  Once 
a  set  of  parameters  is  selected,  the  corresponding  distribution  is  normalized  to  i.e.  unity 
as  g(z)  =  j+^g(zjdz  •  The  bunch  charge  is  then  set  to  Q  =  1  nC  to  yield  the  current  profile 

I(z)  =  Qcg(z)  where  c  is  the  velocity  of  light. 

3.1.1  Simulation  Methods 


In  order  to  explore  the  performances  of  the  current  profiles  in  beam-driven  dielectric- 
wakefield  acceleration,  we  considerer  a  diamond  DLW  with  parameters  listed  in  Tab.  2. 
Given  the  current  profile  /(z),  the  axial  electric  field  is  computed  as  the  integral 


E(z)  =  J2  f  I(^z~  z)Wn(z)dz 

n 


(26) 
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where  Wn(z )  are  the  Green’s  functions  associated  to  the  nth  mode;  see  Ref.  .  For  our 
calculation  with  limit  the  summation  to  n  —  4  modes.  Once  the  axial  field  is  obtained 
the  the  decelerating  field  is  computed  as  EL  =  ma x[E(L)]  for  z  within  the  bunch  and 
E+  =  min[E(2:)]  for  z  behind  the  bunch. 

Table  2:  Parameters  associated  to  the  dielectric  structure  used  in  the  wakefield  simulations. 


Distribution 

Parameters 

Units 

inner  radius  a 

165 

jim 

outer  radius  b 

195 

fim 

relative  permittivity  er 

5.7 

- 

fundamental  frequency  /o 

0.83 

THz 

These  semi-analytical  simulations  were  imbedded  in  a  genetic  optimizer  .  For  each 
current  profile,  the  associated  parameters  were  varied.  The  two  goals  of  the  optimizer  are 
to  find  parameters  that  maximize  E+  and  7Z. 

3.1.2  Results 

Each  shape  was  was  optimized  over  the  parameters  listed  above  respectively.  Bunch  shapes 
with  a  small  number  of  parameters  (e.g.  gaussian  and  ramped  bunch)  converged  more 
quickly  than  more  complicated  shapes  (e.g  double  triangle,  skewed  gaussian,  “realistic” 
parameterization).  On  average,  approximately  15,000  runs  were  done  per  shape.  The  best 
achieved  values  for  1Z  and  E+  are  summarized  in  Fig.  10  (right). 

We  clearly  see  a  trade-off  between  1Z  and  E+:  as  expected,  current  profiles  resulting  in 
large  7Z  are  restricted  to  smaller  values  of  E+  and  vice  versa.  The  data  presented  in  Fig.  10 
was  generated  using  a  1-nC  electron  bunch;  increasing  the  charge  would  result  in  higher  E+ 
without  affecting  1Z.  Interestingly,  none  of  the  asymmetric  shapes  investigated  stand  out 
as  a  best  candidate.  A  linear  regression  of  the  best  cases  provides  an  empirical  limit  for  the 
maximum  value  of  1Z\  ma x.[JZ\  ~  400  x  E^0'8.  For  the  chosen  structure  parameters,  value 
of  E  >  10  are  achieved  with  E+  «  100  MV/m.  Again,  increasing  the  charge  to,  e.g.,  5-nC 
would  result  in  0.5-GV/m  field  or  alternatively  could  enable  reaching  higher  transformer 
ratio  values  for  100-MV/m  fields  as  done  in  Ref.  .  In  DLW-based  acceleration,  larger  E+ 
values  could  also  be  reached  by  reducing  the  aperture  of  the  structure. 

For  the  chosen  structure  parameters,  value  of  1Z  >  10  are  achieved  with  E+  in  the 
vicinity  of  100  MV/m  for  1-nC.  Increasing  the  charge  to,  e.g.,  5  nC  would  results  in  0.5 
GV/m  field  or  alternatively  could  enable  reaching  higher  transformer  ratio  value  for  100- 
MV/m  fields. 

In  Figures.  11  and  12  we  present  an  example  of  optimized  Fourier  and  skew-Gaussian 
distributions,  both  have  1Z  >  6  and  E+  >  150  MV/m.  In  Fig.  12  a  skewed  gaussian  excites 
a  wakefield  with  a  nice  high  1Z.  An  important  aspect  of  the  skew-Gaussian  distribution  on 
this  example  is  its  ability  to  provide  a  slightly  flatted  accelerating  field  over  the  prospective 
location  of  a  witness  bunch.  This  latter  feature,  also  observed  for  the  double-triangle 
distribution,  would  reduce  the  energy  spread  imparted  on  the  witness  bunch. 
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Figure  10:  Shapes  considered  in  our  study:  Gaussian  (blue),  linearly-ramped  (green), 
exponential  (cyan),  Fourier  (red),  double  triangle  (yellow),  and  skewed  gaussian  (magenta) 
(left)  and  trade-off  curve  between  1Z  and  E+  for  the  different  shapes  shown  in  the  left  plot 
with  same  color  coding  (right).  The  axis  of  the  left  plot  are  normalized  to  unity  for  clarity, 
e.g.,  the  traces  shown  correspond  to  #[z/max(z)]/max[g(z)]. 


Figure  11:  Fourier  distribution  (green  trace)  and  associated  wakefield  (blue  trace). 

Finally,  we  present  an  example  of  wakefield  generated  by  the  “realistic”  bunch  in  Fig.  12. 
Such  a  bunch  shape  is  achievable  using  a  dual-frequency  linear  accelerator  which  will 
be  used  in  a  forthcoming  experiment  to  demonstrate  beam-driven  acceleration  with  a 
enhanced  transformer  ratio  . 
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200 


Figure  12:  Skewed-Gaussian  (left)  and  “realistic”  (experimentally  measured)  distributions 
(right)  (green  trace)  and  associated  wakefields  (blue  trace). 

3.2  Optimum  electron-bunch  shape 

In  beam  driven  techniques,  a  “drive”  bunch  passes  through  a  high-impedance  medium  and 
experiences  a  decelerating  field.  The  resulting  energy  loss  can  be  transferred  to  a  properly 
delayed  “witness”  bunch  trailing  the  drive  bunch  in  the  medium.  The  transformer  ratio, 
defined  as  7 Z  =  E+/E-  [where  E+  (reps.  E-)  is  the  maximum  accelerating  (reps,  mini¬ 
mum  decelerating)  field  behind  (reps,  within)  the  drive  bunch]  is  an  important  figure  of 
merit  as  it  describes  the  energy  transfer  within  the  two  bunches.  For  symmetric  bunches, 
1Z  <  2.  Going  beyond  1Z  =  2  requires  asymmetric  bunches  and  can  be  maximized  by 
generating  a  flat  decelerating  field  over  the  drive  bunch;  however,  increasing  7Z  generally 
reduces  E+  and  the  relationship  between  them  is  left  to  a  compromise  . 

In  practice  however,  the  acceleration  gradients  and  acceleration  frequencies  depend  on 
the  interaction  medium,  which  for  e.g.  dielectric  wakefield  acceleration,  depend  on  the  di¬ 
mensions  of  the  structure.  Therefore  altogether,  the  structure  dimensions  are  also  limited 
by  the  beam  properties  such  as  the  beam  energy  and  beta  functions.  At  low  energy  for 
example,  due  to  the  difficulty  of  fitting  into  relatively  long  structures,  it  is  more  valuable 
to  excite  high-gradient  wakes  with  shorter  bunches  with  poor  transformer  ratios.  On  the 
other  hand,  high  energy  beams  with  relatively  smaller  betatron  functions  can  fit  into  much 
longer,  and  smaller  structures.  In  this  perspective  it  is  more  useful  to  drive  wakes  efficiently 
with  high  transformer  ratios. 

To  date,  several  current  profiles  capable  of  generating  optimal  transformer  ratios  have 
been  proposed.  These  include  linearly  ramped  profiles  combined  with  a  door-step  or  ex¬ 
ponential  initial  distribution  ,  and  more  recently  the  double-triangle  current  profiles  .  In 
general,  these  shapes  present  discontinuities  and  rely  on  complicated  beam  manipulation 
techniques  which  have  limited  applications  . 

During  our  research  we  introduce  a  new  class  of  smooth  current  profiles  which  also 
lead  to  constant  decelerating  fields  across  the  drive  bunch  to  lead  to  optimal  transformer 
ratios.  We  also  discuss  a  possible  scheme  for  realizing  these  shapes  employing  a  shaped 
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photocathode-laser  pulse. 

For  simplicity  we  consider  a  wakefield  structure  that  supports  an  axial  wakefield  de¬ 
scribed  by  the  Green’s  function 


G(z)  =  2kcos(uoz/c),  (27) 

where  k,  is  the  loss  factor  and  c^o  =  2nc/X  with  A  being  the  wavelength  of  the  considered 
mode.  Here  £  is  the  distance  behind  the  source  particle  responsible  for  the  wakefield.  We  do 
not  specialize  to  any  particular  wakefield  mechanism  and  depending  on  the  structure  used 
many  modes  might  be  excited  so  that  the  Green’s  function  would  consists  of  a  summation 
over  these  modes. 

3.2.1  Smooth  current  profiles  for  enhanced  transformer  ratios 

Optimal- 77.  shapes  are  useful  for  generating  flat  decelerating  fields  over  the  bunch;  this 
inherently  reduces  the  energy  spread  incurred  over  the  bunch.  Additionally,  these  shapes 
lead  to  the  most  efficient  energy  transfer  between  a  drive  and  witness  bunch  which  could 
lead  to  overall  higher  witness  bunch  energies. 

Based  on  the  work  presented  in  Ref.  ,  we  consider  a  bunch-current  profile  described  as 
a  piecewise  function  I(t).  We  take  I(t)  to  be  piecewise  on  two  intervals  [0,  r]  and  [r,  T]  and 
zero  elsewhere.  Also  we  constrain  our  search  to  functions  such  that  I{t)  and  I(t)  =  dl/dt 
are  continuous  at  t  —  r. 

As  an  example  with  consider  a  modulated  linear  ramp  function  over  [0,  r]  of  the  form 

s(t)  =  at  +  6sin(o;t),  (28) 


where  a  and  b  are  constant,  and  c^o  is  the  modulation  frequency.  We  correspondingly 
introduce  the  current  profile  as  the  piecewise  ) continuous)  function 


m  =  j0 


s{t) 

s'(r)t  —  s'(t)t  +  s(t) 

0 


if  0  <  t  <  r, 
if  t  <  t  <  T, 
elsewhere. 


(29) 


where  Iq  is  the  peak  current.  An  example  of  current  shape  is  depicted  in  Fig.  13  (a). 


We  further  specialize,  for  simplicity,  to  the  case  where  r  —  2u7r/cu  (where  n  G  N)  and 
—  uo/q  (where  cjq  is  defined  in  Eq.  27  and  m  G  N)  )  The  decelerating  field  (H-)  over 
the  drive  bunch  is  found  to  be 


~U2$-1)  [}aq2  -  (bu  +  a)]  COS  Y  +  bu  cos Ldt  -  a(q2  -  1)) 


V~(t)  =  I0{ 


if  0  <  t  <  t 


^2  (g2 _ 1}  (  buq2  cos 


ujt—27rn 


o 


+  [ aq 2  —  (bco  +  a)]  cos  y  +  (bu  +  a)(l  —  q2)^j  if  0  <  t  <  r 

elsewhere. 

(30) 
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Figure  13:  Example  (a)  of  current  profile  (blue  trace)  described  by  Eq.  29  (a)  and  s(t) 
function  (red  trace).  The  corresponding  induced  voltage  (b)  within  (blue  trace)  and  behind 
(red  trace)  the  bunch.  For  these  data  Iq  =  1,  a  =  1,  b  =  3a/5,  n  =  1,  m  =  5  and  q  —  2. 
(the  current  and  voltage  units  are  arbitrary). 


The  above  equation  shows  that  when  ^  —  \  (where  l  G  N),  the  decelerating  field  is  constant 
for  t  G  [t,  T)  and  can  be  written 


Vte[r,T\^)  =  ~2 -Q  ([Mt”1)  M2  -  1  )  +  a{q2  -  l)]cos^  +  (bcu  +  a)(l  -  q2)^j  , 

1  V  (31) 

so  that  the  time  dependence  of  the  decelerating  potential  can  be  cancelled  by  choosing  the 
ramp  parameters  a  and  b  to  satisfy  the  equation 


a(  1  ~  g2) 

u[(-l)Tq2-l] 


(32) 


Specializing  to  the  case  where  2 n/q  is  an  odd  number  gives  b  =  a(q 2  —  1)/\uj(1  +  q2)\  and 
the  decelerating  voltage  across  the  bunch  takes  the  form 


v-(t)  =  i0{ 


aq2Io 

oj2(l+q2) 
2  alog4 

0J2(l+q2) 

0 


1  cos  y  +  a  cos  wt  —  (1  +  q* 


if  0  <  t  <  r, 
if  r  <  t  <  T, 
elsewhere. 


(33) 


The  latter  case  is  illustrated  in  Fig.  13  (b)  with  l  =  n/q  =  1/2.  Likewise  the  induced 
oscillating  voltage  behind  the  bunch  can  be  obtained  from 

V+(t)  =  2(1+  2)  —  2m)  —  1]  S^n  —  nm'j  +  2q2  cos  —  7rm^  (34) 

The  second  term  can  be  dropped  to  good  approximation  to  find  a  lower  estimate  on 
the  transformer  ratio,  without  loss  of  generality 

n  =  (m + 7r’  (35) 
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where  m  is  the  number  of  the  fundamental-mode  wavelengths  within  the  total  bunch  length. 
We  note  that  for  short  bunch  durations  (i.e.  as  needed  to  produce  high  accelerating  fields), 
the  latter  transformer  ratio  is  higher  than  other  proposed  distributions. 

3.2.2  Beam  shaping  using  a  tailored  photocathode  laser 

To  date,  longitudinal  shaping  for  enhanced  transformer  ratios  has  only  been  demonstrated 
using  a  mask  in  a  dispersive  section.  Another  proposed  technique  relies  on  using  a  mask 
to  shape  a  beam  transversely,  and  use  an  emittance  exchanger  to  rotate  the  shape  longi¬ 
tudinally.  In  the  former  case,  the  scheme  is  limited  to  specific  shapes  in  the  mask,  and 
inherits  beam  losses.  Additionally  shapes  which  require  sharp  edges  (e.g.  double-triangle) 
are  limited  by  the  dispersion  and  emittance  of  the  beam  at  the  mask.  The  latter  scheme 
requires  a  deflection  cavity  in  a  dispersive  section  which  complicates  the  overall  procedure; 
moreover,  this  scheme  also  inherits  beam  losses  from  a  mask  upstream,  and  is  also  limited 
by  beam  emittances. 

An  alternative  approach  is  to  longitudinally  tailor  the  laser  profile  in  a  high-gradient 
photoinjector.  In  such  a  scenario,  electrons  emitted  from  the  photocathode  will  have  an 
initial  longitudinal  distribution  similar  to  the  injected  laser  pulse.  We  consider  an  example 
of  an  S-Band  gun  operating  at  140  MV/m  to  reduce  the  space-charge  effects.  Generally 
higher  accelerating  fields  in  a  gun  lead  to  higher  charge  densities  which  could  drive  larger 
gradient  wakefields  inside  DLWs.  We  use  the  particle-tracking  program  astra  which  takes 
into  account  space  charge  effects  using  a  cylindrical-symmetric  algorithm  .  We  model  the 
1-nC  electron  bunch  considered  in  our  studies  with  50,000  macroparticles.  Additionally 
the  laser  is  chosen  to  have  a  transverse  rms  spot  size  of  ac  =  0.8  mm  and  rms  duration  of 
(Jt  —  1  ps.  A  solenoidal  lens  is  placed  downstream  of  the  gun  to  maintain  the  transverse 
size  below  1  cm  over  a  drift  length  of  1  m. 

In  general,  photoemission  is  a  intricate  process  which  depends  on  many  parameters 
including  the  amplitude  and  phase  (with  respect  to  the  laser)  of  the  applied  accelerating 
field  on  the  photocathode  surface,  the  bunch  charge,  and  the  cathode  material  properties. 
It  is  challenging  to  analytically  find  an  optimum  laser  shape  that  provide  a  given  current 
distribution  at  a  given  location.  Work  toward  the  production  of  ramped  current  profile 
was  investigated  in  the  context  of  seeded  free-electron  laser  research  and  experimentally 
demonstrated  at  the  FERMI  facility  .  Additionally,  we  demonstrated  that  a  passive  tech¬ 
nique  employing  a  wakefield  structure  combined  with  a  non-ultra-relativistic  bunch  could 
produce  ramped  bunch. 

Therefore  to  explore  the  production  of  tailored  electron  bunches  using  temporally- 
shaped  photocathode  laser  pulses,  we  carry  numerical  simulation  using  the  ASTRA  particle 
tracking  program  .  and  explored  the  outcome  of  various  shapes.  In  this  section  we  report 
on  the  performance  of  the  exponential  distribution  given  by: 

f(t)  =  /0(e^  -  1  )H(6  -  t ),  (36) 

where  fi  is  a  growth  rate  and  6  is  the  ending  time  of  the  pulse  H{t)  is  the  Heaviside  function. 
The  laser  intensity  /o  has  to  be  normalized  to  produce  the  required  bunch  charge  after  the 
photo-emission  process.  We  consider  three  case  of  exponentially-shaped  laser  distributions. 
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All  have  a  total  pulse  duration  of  9  =  5  ps  but  different  steepness  characterized  by  a  growth 
rate  (i  —  3.90  x  105,  6.92  x  105,  1.00  x  106;  see  Fig.  14  (upper-right  plot). 

The  growth  rate  fi  primarily  influences  the  electrostatic  field  inside  the  bunch  which 
affects  the  final  longitudinal  phase  space.  The  case  of  large  growth  rate  yield  the  formation 
of  ellipsoidal  distribution  via  the  blow-out  regime.  In  the  limit  of  small  growth  rates,  the 
laser  profile  become  linearly  ramped  but  the  ramp  quickly  dissipates.  Figure  14  (upper 
right,  and  lower  row  plots)  displays  snapshots  of  the  the  longitudinal  phase  spaces  and 
associated  current  distribution  for  the  three  case  of  growth  rate  mentioned. 


Figure  14:  Initial  photocathode  laser  pulse  temporal  distributions  (top  left)  and  snapshots 
of  the  current  (a)  and  longitudinal  phase  space  (b)  at  25  cm  (red),  50  cm  (green),  100  cm 
(blue)  from  the  photocathode.  The  top-right,  lower-left  and  lower-right  plots  correspond 
respectively  to  the  initial  laser  distributions  appearing  as  a  red,  green,  and  blue  trace  (with 
increasing  steepness)  in  the  top-left  plot. 


Finally,  we  quantify  the  capability  of  the  formed  electron  bunches  to  serve  as  a  wake- 
field  drive  by  considering  their  applications  to  drive  wakefield  in  a  dielectric  structure.  We 
therefore  convolve  the  shapes  with  the  one-dimensional  Green’s  function  associated  to  a 
cylindrically-symmetric  dielectric-line  waveguide  (DLW)  .  We  choose  a  DLW  with  inner 
radius  a  —  165  /im,  outer  radius  b  =  195  /im,  and  relative  dielectric  permittivity  e  =  5.7; 
here  we  note  that  in  order  to  fit  the  beam  into  a  structure  with  such  dimensions  would 
require  a  significant  increase  in  beam  energy  which  could  also  “freeze”  the  evolution  of  the 
beam.  For  the  three  cases,  we  record  E+  and  1Z  along  the  drift  following  the  photocath¬ 
ode,  the  results  are  shown  in  Fig.  15.  The  results  indicate  that  both  large  E+  and  1Z  are 
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Figure  15:  The  calculated  E+  (blue  trace)  and  1Z  (red  trace)  as  a  function  of  drift  length 
from  the  photocathode  in  a  DLW  with  dimensions  described  above.  The  top  left,  right, 
and  and  bottom  left  plots  correspond  to  initial  exponential  distributions  red,  green,  blue 
shown  in  Fig.  14  (top  left)  respectively. 


achievable  for  the  different  laser  distribution  along  the  beamline.  The  variation  parameter 
/i  adds  some  flexibility  between  large  E+  and  R  and  generally  larger  /i  lead  to  larger  E+. 
Although  in  our  examples,  the  bunch  shape  is  still  evolving,  it  could  in  principle  be  possible 
to  add  a  booster  cavity  and  optimize  the  setup  to  ensure  the  optimum  shape  is  formed 
downstream  of  the  booster  cavity  at  a  sufficiently  high  energy  to  freeze  the  longitudinal 
motion  within  the  bunch. 
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3.3  Experimental  realization  of  a  linearly- ramped  bunch  with 

A  MULTIFREQUENCY  LINAC 

We  first  elaborate  the  proposed  method  using  a  1D-1V  single-particle  model  of  the  LPS 
dynamics  and  take  an  electron  with  coordinates  (z,  5)  where  z  refers  to  the  longitudi¬ 
nal  position  of  the  electron  with  respect  to  the  bunch  barycenter  (in  our  convention 
z  >  0  corresponds  to  the  head  of  the  bunch)  and  5  =  p/(p)  —  1  is  the  fractional  mo¬ 
mentum  error  (p  is  the  electron’s  momentum  and  (p)  the  average  momentum  of  the 
bunch).  Considering  a  photo-emission  electron  source,  the  LPS  coordinates  downstream 
are  (zo,5o  —  cloZq  +  ^o+j  +  ®(zo))  where  ao,  and  &o  are  constants  that  depend  on  the 
bunch  charge  and  operating  parameters  of  the  electron  source.  For  sake  of  simplicity 
we  limit  our  model  to  second  order  in  zq  and  Next,  we  examine  the  acceleration 
through  a  linac  operating  at  the  frequencies  fi  and  fn  =  nfi  with  total  accelerating  volt¬ 
age  V(z )  =  V±  cos(kiz  +  Lp\)  +  Vn  cos (knz  +  cpn)  where  Vyn  and  cpijU  are  respectively  the 
accelerating  voltages  and  operating  phases  of  the  two  linac  sections,  and  k\^n  =  27r/i>n/c. 
In  our  convention,  when  the  phases  between  the  linac  sections  and  the  electron  bunch 
are  <piin  —  0  the  bunch  energy  gain  is  maximum  (this  is  refer  to  as  on-crest  opera¬ 
tion).  Under  the  assumption  kiiUzo  <C  1  and  neglecting  non-relativistic  effects,  the  elec¬ 
tron’s  LPS  coordinate  downstream  of  the  linac  are  (zi  =  zo,5i  =  aizo  +  b[Z q)  where 
ai  =  a0  -  e{k\V\  sin^i  +  knVn  sin  tpn )/Eu  bt  =  b0  -  e(k\V\  cos^i  +  k%Vn  cos  tpn) / (2Et) 
with  e  being  the  electronic  charge  and  Ei  the  beam’s  average  energy  downstream  of  the 
linac.  Finally,  we  study  the  passage  of  the  bunch  through  an  achromatic  current-enhancing 
dispersive  section  [henceforth  referred  to  as  “bunch  compressor”  (BC)].  The  LPS  dynam¬ 
ics  through  a  BC  is  approximated  by  the  transformation  Zf  —  R^Si  +  T^Sf  where 
(also  referred  to  as  longitudinal  dispersion),  and  T$qq  are  the  coefficients  of  the  Taylor 
expansion  of  the  transfer  map  (zf\Si)  of  the  BC.  Therefore  the  final  position  is  given 
as  function  of  the  initial  coordinates  following  Zf  =  afZo  +  bfZ^  with  aj  =  1  +  aiR§§ 
and  bf  =  biR$Q  +  afT^QQ.  Taking  the  initial  current  to  follow  the  Gaussian  distribu¬ 
tion  Io(zo)  =  Iq  exp[— Zq/{2g2z  0)]  (where  Iq  is  the  initial  peak  current),  and  invoking  the 
charge  conservation  If(zf)dz  =  Io(zo)dzo  gives  the  final  current  distribution  IJ(zf)  = 

/0/A1/2(z/)exp[— (af  +  A1/2(zf))2/(8b2fal0)}@[A(zf)]  where  A (zf)  =  aj  +  AbfZf  and  ©() 
is  the  Heaviside  function.  The  latter  current  distribution  does  not  include  the  effect  of  the 
initial  uncorrelated  fractional  momentum  spread  a%0.  The  final  current,  taking  into  ac¬ 
count  cr^0,  is  given  by  the  convolution  If(zf)  =  f  dzflj(zf)  exp  [—(zf  —  i/)2/(2cr^)]  where 
<ju  =  i?56^o*  currenf  shape  is  controlled  via  aj  and  bf  and  can  be  tailored  to 

follow  a  linear  ramp  as  demonstrated  in  Fig.  16. 

The  experiment  described  in  this  Letter  was  performed  at  the  Free-electron  LASer 
in  Hamburg  (FLASH)  facility  .  In  the  FLASH  accelerator,  diagrammed  in  Fig.  17,  the 
electron  bunches  are  generated  via  photoemission  from  a  cesium  telluride  photocathode 
located  on  the  back  plate  of  a  1+1/2  cell  normal-conducting  rf  cavity  operating  at  1.3 
GHz  on  the  TMqio  7r-mode  (rf  gun).  The  bunch  is  then  accelerated  in  a  1.3-GHz  and  3.9- 
GHz  superconducting  accelerating  modules  (respectively  ACC1  and  ACC39)  before  passing 
through  a  bunch  compressor  (BC1).  The  ACC39  3rd-harmonic  module  was  installed  to 
nominally  correct  for  nonlinear  distortions  in  the  LPS  and  enhance  the  final  peak  current 
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Figure  16:  Analytically-computed  current  profiles  for  several  values  of  bf  for  fixed  af  =  2.5 
(a)  and  for  several  values  of  af  with  bf  =  0.7  (b).  The  numbers  in  (a)  [resp.  (b)]  are  the 
values  of  bf  [resp.  aj];  for  all  the  cases  au  =  0.05. 


of  the  electron  bunch  .  Downstream  of  BC1,  the  bunch  is  accelerated  and  can  be  further 
compressed  in  BC2.  A  last  acceleration  stage  (ACC4/5/6/7)  brings  the  beam  to  its  final 
energy  (maximum  of  ~  1.2  GeV).  The  beam’s  direction  is  then  horizontally  translated 
using  a  dispersionless  section  referred  to  as  dogleg  beamline  (DLB).  Nominally,  the  beam 
is  sent  to  a  string  of  undulators  to  produce  ultraviolet  light  via  the  self-amplified  stimulated 
emission  free-electron  laser  (FEL)  process.  For  our  experiment,  the  bunches  were  instead 
vertically  sheared  by  a  2.856-GHz  transverse  deflecting  structure  (TDS)  operating  on  the 
TMno-like  mode  and  horizontally  bent  by  a  downstream  spectrometer  .  Consequently 
the  transverse  density  measured  on  the  downstream  Cerium-doped  Yttrium  Aluminum 
Garnet  (Ce:YAG)  scintillating  screen  is  representative  of  the  LPS  density  distribution. 
The  horizontal  and  vertical  coordinates  at  the  Ce:YAG  screen  are  respectively  xs  ~  77 
where  77  ~  0.75  m  is  the  horizontal  dispersion  function,  and  ys  ~  KZp  where  ft  ~  20  is  the 
vertical  shearing  factor  and  (zf,5f)  refers  to  the  LPS  coordinate  upstream  of  the  TDS. 
The  exact  values  of  77  and  ft  are  experimentally  determined  via  a  beam-based  calibration 
procedure. 

The  accelerator  parameters  settings  are  gathered  in  Tab.  3.  The  nominal  settings 

(2) 

of  BC2  were  altered  to  reduce  its  longitudinal  dispersion  R^6J  and  the  ACC2/3  and 
ACC4/5/6/7  accelerating  modules  were  operated  on  crest.  Such  settings  insure  that  the 
BC2  and  the  DBL  sections  do  not  significantly  affect  the  LPS  beam  dynamics.  Therefore 
the  measured  current  profile  is  representative  of  the  profile  downstream  of  BC1. 

In  order  to  validate  the  simple  analytical  model  described  above,  numerical  simulations 
of  the  LPS  beam  dynamics  were  carried  using  a  multi-particle  model.  The  simulations  also 
enable  the  investigation  of  possible  detrimental  effects  resulting  from  collective  effects  such 
as  longitudinal  space  charge  (LSC)  and  beam  self  interaction  via  coherent  synchrotron 
radiation  (CSR)  .  In  these  simulations,  the  beam  dynamics  in  the  rf-gun  was  modeled  with 
the  particle-in-cell  program  ASTRA  and  the  obtained  distribution  was  subsequently  tracked 
in  the  accelerating  modules  using  a  ID- IV  program  that  incorporates  a  one-dimensional 
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Table  3:  Settings  of  accelerator  subsystems  relevant  to  the  LPS  dynamics  used  in  the 
experiment  and  simulations. 


parameter 

symbol 

value 

unit 

AC  Cl  voltage 

u 

[140-157] 

MV 

ACC1  phase 

¥>i 

[-10,10] 

deg 

ACC39  voltage 

[13,21] 

MV 

ACC39  phase 

¥>3 

[160-180] 

deg 

ACC2/3  voltage 

U,2-3 

311 

MV 

ACC2/3  phase 

¥>1,2-3 

0 

deg 

ACC4/5/6/7  voltage 

U,4— 7 

233.9 

MV 

ACC4/5/6/7  phase 

¥>1,4—7 

0 

deg 

BC1  longitudinal  dispersion 

p(l) 

-^56 

~  170 

mm 

BC2  longitudinal  dispersion 

p(2) 

~  15 

mm 

Single-bunch  charge 

Q 

0.5 

nC 

Bunch  energy 

E 

~  690 

MeV 

rf9Un  Rfil 

ACC1 

_ ]\nv  Ny 


□  c 


X  / 

laser  ACC39 


ACC2/3 


BC2 


LPS  diagnostics^Ce:YAG  screen 
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Figure  17:  Diagram  of  the  FLASH  facility.  Only  components  affecting  the  longitudinal 
phase  space  beam  (LPS)  dynamics  of  the  bunches  are  shown.  The  acronyms  ACC,  BC,  and 
DBL  stand  respectively  for  accelerating  modules,  bunch  compressors,  and  dogleg  beamline 
(the  blue  rectangles  represent  dipole  magnets).  The  transverse  deflecting  structure  (TDS), 
spectrometer  and  Ce:YAG  screen  compose  the  LPS  diagnostics. 


model  of  the  LSC.  The  program  CSRTRACK  ,  which  self-consistently  simulates  CSR  effects, 
was  used  to  model  the  beam  dynamics  in  the  BC1,  and  BC2  sections.  An  example  of 
simulated  LPS  distributions  and  associated  current  profiles  computed  for  different  settings 
of  ACC1  and  ACC39  parameters  appear  in  Fig.  18.  The  results  indicate  that  the  produc¬ 
tion  ramped  bunches  is  possible  despite  the  intricate  LPS  structures  developing  due  to  the 
collective  effects  and  higher-order  nonlinear  effects  not  included  in  our  analytical  model. 
The  simulations  also  confirm  that  the  current  profile  upstream  of  the  TDS  (as  measured 
by  the  LPS  diagnostics)  is  representative  of  the  one  downstream  of  BC1. 

Figure  19  displays  examples  of  measured  LPS  distributions  with  associated  current 
profiles  obtained  for  different  settings  of  ACC1  and  ACC39.  As  predicted,  the  observed 
current  profiles  are  asymmetric  and  can  be  tailored  to  be  ramped  with  the  head  of  the 
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Figure  18:  Simulated  LPS  distribution  [(a)  and  (b)]  with  associated  current  profile  down¬ 
stream  of  BC1  (solid  blue  trace)  and  DBL  (dash  red  trace)  [(c)  and  (d)].  The  set  of  plots 
[(a),  (c)]  and  [(b),  (d)]  correspond  to  different  (V 1,3,  ^1,3)  settings. 


bunch  (z  >  0)  having  less  charge  than  the  tail;  see  Fig.  19  (b-d).  The  latter  feature  is  in 
contrast  with  the  nominal  compression  case  at  FLASH  where  the  LPS  distortion  usually 
results  in  a  low-charge  trailing  population  as  seen  in  Fig.  19  (a). 

We  now  quantify  the  performance  of  the  produced  current  profiles  to  enhance  beam- 
driven  acceleration  techniques  by  considering  a  drive  bunch  injected  in  a  cylindrical- 
symmetric  dielectric-lined  waveguide  (DLW)  .  The  DLW  consists  of  a  hollow  dielectric 
cylinder  with  inner  and  outer  radii  a  and  b.  The  cylinder  is  taken  to  be  diamond  (relative 
electric  permittivity  er  =  5.7);  and  its  outer  surface  is  contacted  with  a  perfect  conduc¬ 
tor;  see  Fig.  20  (a).  The  measured  current  profiles  are  numerically  convolved  with  the 
Green’s  function  associated  to  the  monopole  mode  to  yield  the  axial  electric  field  .  These 
semi-analytical  calculations  were  benchmarked  against  finite-difference  time-domain  elec¬ 
tromagnetic  simulations  executed  with  VORPAL  .  The  transformer  ratio  is  numerically 
inferred  as  7Z  =  \E+/E-\  where  E_  (resp.  E+)  is  the  decelerating  (resp.  accelerating)  ax¬ 
ial  electric  field  within  (resp.  behind)  the  electron  bunch;  see  Fig.  20  (b).  The  achieved  1Z 
and  E+  values  as  the  structure  geometry  is  varied  are  shown  in  Fig.  21.  As  a  G  [20,  300]  /im 
and  b  G  a+  [20,  300]  /im  are  varied  the  wavelengths  of  the  excited  wakefield  modes  change. 
The  simulations  show  that  profiles  (b)  and  (c)  of  Fig.  19  can  yield  values  of  1Z  >  2.  A 
possible  configuration  with  (a,  b)  —  (20,  60)  /im,  results  in  1Z  ~  5.8  with  E+  cz  0.75  GV/m; 
see  corresponding  wake  in  Fig.  20  (b).  Such  high-field  with  transformer  ratio  significantly 
higher  than  2  and  driven  by  bunches  produced  in  a  superconducting  linac  could  pave  the 
way  toward  compact  high-repetition-rate  short-wavelength  FELs  . 

Finally,  the  proposed  technique  could  be  adapted  to  non-ultrarelativistic  energies  using 
a  two-  (or  multi-)  frequency  version  of  the  velocity-bunching  scheme  .  Such  an  implemen¬ 
tation  would  circumvent  the  use  of  a  BC  and  would  therefore  be  immune  to  CSR  effects. 
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Figure  19:  Snapshots  of  the  measured  longitudinal  phase  spaces  (left  column)  and  as¬ 
sociated  current  profiles  (right  column)  for  different  settings  of  the  ACC1  and  ACC39 
accelerating  modules.  The  values  (Vi, ipi\ V3, <^3)  [in  (MV,°;  MV,0)]  are:  (150.5,  6.1;  20.7, 
3.8),  (156.7,  3.8;  20.8,  168.2),  (155.6,  3.6;  20.6,  166.7),  and  (156.8,  4.3;  20.7,  167.7)  for 
respectively  case  (a),  (b),  (c),  and  (d). 
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Figure  20:  Cylindrical-symmetric  dielectric-loaded  waveguide  considered  (a)  and  axial 
wakefield  produced  by  the  current  profile  shown  in  Fig.  19  (c)  for  (a,  6)=(20,60)  /tin. 


4  Experimental  setup  for  slab-dielectric  structure 

TEST  AT  FERMILAB 
4.1  EXPERIMENTAL  SETUP 

The  beamline  configuration  for  our  DWFA  experiment  is  diagrammed  in  Fig.  22.  The 
beamline  comprises  a  L-band  RF  gun  followed  by  two  SCRF  accelerating  cavities  (CAV1 
and  2).  The  RF  gun  is  nested  in  a  pair  of  solenoidal  lenses  that  can  be  used  to  produce 
beam  with  large  angular  momentum.  Such  a  beam  can  be  decoupled  by  a  set  of  three 
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transformer  ratio  a  (|Lim) 

Figure  21:  Simulated  A  ratio  versus  peak  accelerating  field  (a)  for  the  four  measured 
current  profiles  (displayed  as  different  colors  with  label  corresponding  to  cases  shown  in 
Fig.  19).  Transformer  ratio  (false  color  map)  as  a  function  of  the  DLW  inner  radius  a 
and  dielectric  layer  thickness  b  —  a  with  corresponding  \E+\  shown  as  isoclines  with  values 
quoted  in  MV/m  for  case  (c)  of  Fig.  19. 


skew  quadrupole  magnets  downstream  of  CAV2  to  produce  flat  beams  —  beams  with  high 
transverse  emittance  ratio  —  that  can  then  be  compressed  using  a  magnetic  chicane  (BC1). 
The  skew  quadrupole  magnets  insertion  is  referred  to  as  round-to-flat-beam  transformer 
(RFBT).  Downstream  of  BC1,  a  triplet  is  used  to  focus  the  beam  inside  the  DLW  structure 
mounted  in  a  2-axis  goniometer;  additionally,  a  linear  stages  will  give  control  over  the 
aperture  of  the  DLW  by  varying  the  gap  between  the  two  slabs.  The  beam  is  finally  drifted 
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Figure  22:  Layout  of  the  ASTA  photoinjector.  LI  and  L2  are  solenoids,  and  CAV1  and 
2  SCRF  cavities.  The  green  and  re  rectangles  respectively  represent  normal  and  skewed 
quadrupole  magnets.  X124  is  Ce:YAG  a  screen  located  in  the  dispersive  section. 

to  the  vertical  spectrometer  and  directed  to  a  Ce:YAG  screen  (X124)  with  full  vertical  size 
of  38  mm  and  located  at  a  dispersion  of  r\y  —  0.44  m  permitting  the  measurement  of  beam 
spectrum  of  5p/p  ~  9  %  relative  momentum  spread.  The  high-resolution  CCD  (Prosilica 
GC  2450)  could  in  principle  enable  resolution  below  Sp/p  ~  10-4  for  an  ideal  zero-emittance 
beam.  The  zero-charge  betatron  functions  at  X124  viewer  are  shown  in  Fig.  23  as  a  function 
of  betatron  functions  obtained  at  the  waist  (ax  =  ay  =  0)  in  the  DLW  structure.  The 
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focusing  between  the  DLW  structure  and  X124  screen  is  solely  achieved  by  the  dipole  (no 
quadrupole  magnets  are  presently  installed  in  this  section).  For  a  vertical  beta  function 


P*ip  (m)  ftp  (m) 


Figure  23:  The  zero-charge  betatron  functions  (left)  and  the  transverse  RMS  beam  sizes 
(right)  along  the  ASTA  accelerator  with  flat  beams. 

of  /3*  ~  2  m  at  the  center  of  the  DLW,  the  resulting  f3  function  at  X124  is  /3;f 124  ~  1  m 
limiting  the  energy  resolution  of  the  spectrometer  to  1.8  x  10-4  (for  a  geometric  emittance 
of  ~  6.3  nm  (corresponding  to  0.5  /am  normalized  with  a  Lorentz  factor  7  ~  80). 

Finally,  a  diagnostics  station  located  downstream  of  the  vertical  spectrometer  will  en¬ 
able  the  detection  and  autocorrelation  of  THz  radiation  generated  by  the  bunch  passing 
through  the  DLW  structure.  The  autocorrelation  is  performed  with  a  Michelson  interfer¬ 
ometer  employing  a  He-cooled  InSb  bolometer;  see  details  in  Ref.  . 

Our  experiment  relies  on  the  production  of  a  flat  beam,  i.e.,  a  beam  with  large  trans¬ 
verse  emittance  ratio  .  In  our  setup  we  produce  flat  beams  with  a  low  vertical  emittance  to 
mitigate  horizontal-emittance-dilution  effects  arising  in  BC1  via  space  charge  and  coherent 
synchrotron  radiation.  Another  benefit  of  this  configuration  is  the  low  betatron  contribu¬ 
tion  to  the  beam  size  at  X124  given  a  vertical  normalized  emittance  as  low  as  ey  ~  0.3  fi m. 
An  important  challenge  to  overcome  is  the  formation  of  compressed  flat  beam  as  described 
elsewhere  . 

4.2  Start-to-end  simulations 

The  start-to-end  simulations  detailed  below  were  performed  using  particle-in-cell  beam- 
dynamics  program  including  ASTRA  and  IMPACT-T  .  The  distribution  downstream  of  the 
compressor  was  then  matched  to  a  waist  at  the  DLW  structure  location  with  ELEGANT  .  To 
model  the  beam  self-interaction  with  its  wakefield  in  the  DLW,  we  use  a  modified  version  of 
IMPACT-T  described  in  Ref.  .  The  dielectric-wakefield  model  is  based  on  a  3-D  Green’s  func¬ 
tion  approach  and  was  successfully  benchmarked  against  particle  in  cell  finite-difference 
time-domain  simulations  .  We  consider  a  DLW  structure  composed  of  two  parallel  dielec¬ 
tric  slabs.  The  separation  between  the  inner  surface  and  outer  (metallized)  surfaces  is 
respectively  2 a  and  2b.  The  dielectric  thickness  is  b  —  a  and  its  relative  permittivity  is 
taken  to  be  er  —  5.7  to  correspond  to  diamond. 

4.2.1  Case  of  single-mode  DLW  structures 

Single-mode  structures  have  the  advantage  to  produce  sinusoidal  fields  with  known  wave¬ 
lengths.  However  since  the  beam’s  energy  couples  to  a  single  mode,  the  resulting  accel¬ 
erating  fields  are  generally  smaller  than  the  accelerating  gradients  achieved  in  multimode 
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structures. 


Figure  24:  Longitudinal  phase  spaces  (a,c)  and  associated  transverse  distributions  at  X124 
(b,d)  with  (a,b)  and  without  (c,d)  a  DLW  structure.  The  DLW  structure  used  for  (c,d) 
has  parameters  a=100  /im,  b  =  120  /im,  and  er  =  5.7 

Due  to  the  limited  measurement  apparatus  in  the  planned  experiment,  we  must  es¬ 
tablish  a  method  of  decoding  the  information  from  the  projected  energy  modulated  beam 
onto  X124.  Figure  24(a,b)  illustrates  the  nominal  longitudinal  phase  space  (LPS)  and  the 
associated  transverse  beam  distribution  at  the  XI 24  location  when  no  DLW  structure  is 
inserted  in  the  beam  path. 

A  complication  arising  from  the  relatively  long  photocathode  laser  pulse  (rms  duration 
at  —  3  ps)regards  the  accumulation  of  a  strong  quadratic  distortion  on  the  LPS  during  ac¬ 
celeration  in  CAV1  and  CAV2  which,  after  bunch  compression  results  in  a  highly  distorted 
distribution.  The  distribution  has  some  benefits  to  the  investigation  of  wakefields  as  the 
charge  concentration  in  the  bunch  head  lead  to  high  peak  current  5  kA)  that  excites 
strong  wakes  while  the  long  trailing  electron  population  samples  this  wake  over  several 
periods.  This  feature  is  clearly  demonstrated  in  Fig.  24(c)  where  the  LPS  immediately 
downstream  of  a  DLW  structure  with  parameters  a=100  /im,  b  =  120  /im,  and  er  —  5.7  is 
shown.  The  resulting  transverse  distribution  on  X124  shows  some  horizontal  bands  that 
correspond  to  the  local  maxima  of  the  observed  energy  modulations  on  the  LPS.  Because 
of  the  large  number  of  modulations,  some  smearing  occurs  at  X124.  In  addition,  we  note 
that  the  Cherenkov  pattern  resulting  from  the  dependence  of  the  accelerating  field  on 
the  transverse  coordinates  can  be  clearly  resolved  at  XI 24  and  could  provide  insightful 
measurements  for  precise  benchmarking  of  the  3-D  model. 

4.2.2  Case  of  multi-mode  DLW  structures 

The  high-peak  current  and  narrow  width  of  the  bunch  head  is  capable  of  exciting  has  a 
spectrum  that  extents  close  to  10  THz.  It  is  therefore  capable  to  excite  higher-order  modes 
possibly  supported  by  the  DLW  structure.  A  simple  way  to  enable  multi-mode  operation 
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of  DLW  structures  consists  of  using  of  thicker  dielectric  thicknesses.  Additionally,  a  larger 
thickness  can  also  lead  to  higher  accelerating  fields  in  the  DLW  as  multiple  modes  can 
constructively  add  up.  The  draw  back  being  the  lack  of  control  on  the  superposition  of  the 
mode  which  can  possibly  lead  to  accelerating  field  region  with  smaller  longitudinal  extend 
compared  to  single-mode  structures.  Such  a  lack  of  control  could  lead  to  the  acceleration 
of  witness  bunches  with  distorted  LPS  or  would  require  very  short  witness  bunches. 

Figure.  25  shows  examples  of  LPS  and  transverse  beam  distribution  simulated  at  X124 
for  two  dielectric  structures  with  inner  radius  a  —  100  /im  and  outer  radius  b  —  150  and 
200  /im.  The  experimental  advantage  for  using  multimode  structures  are  the  lower  number 
of  energy  modulations  which  leads  to  fewer  (and  brighter)  energy  (horizontal)  bands  at 
X124. 


0  0.6  1.2  1.8  -15  -10  -5  0  5  10  15 
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Figure  25:  Longitudinal  phase  spaces  (a,c)  and  associated  transverse  distributions  at  X124 
(b,d)  for  a  dielectric  thickness  of  r  =  50  (a,b)  and  100  jim  (c,d)  DLW  structure.  The  DLW 
structure  other  parameters  are  a  —  100  /im,  b  —  a  +  r,  and  er  —  5.7. 

For  the  case  of  multimode  structures  the  peak  accelerating  field  are  110  and  120  MV/m 
for  respectively  a  dielectric  thickness  of  r  —  50  [Fig.  25(a,b)]  and  100  fim  [Fig.  25(c,d)]  while 
the  single-mode  structure  was  supporting  a  field  amplitude  of  105  MV/m  [Fig.  24(c,d)]. 
We  finally  note  that  using  a  flatter  beam  could  allow  for  smaller  inner-aperture  structure 
which  could  lead  to  larger  accelerating  fields.  Another  path  to  high  axial  field  is  the 
production  of  flat  beams  with  higher  peak  current.  This  latter  possibility  would  rely  on  the 
use  of  an  enhanced  compression  scheme  employing  a  high-harmonic  accelerating  cavity  or 
wakefield  linearizer  to  linearizer  the  bunch’s  LPS  prior  to  the  BC1  chicane.  Ultimately, 
given  limited  vertical  aperture  of  the  X124  screen  (38  mm),  the  measurement  of  larger 
accelerating  gradients  would  rely  on  the  use  of  shorter  DLW  structures  or  lower-charge 
bunches. 


final  report  —  HDTRA-1-10-1-0051 


page  35  of  43 


4.3  Conclusion  &  Outlook 


We  have  demonstrated  that  ASTA  could  support  tests  of  beam-driven  acceleration  based 
on  a  slab-symmetric  DLW  structure  driven  by  flat  beams  leading  to  accelerating  gradients 
larger  than  100  MV/m. 

We  note  that  a  possible  extension  will  be  the  incorporation  of  longitudinally  shaped 
bunches  which  could  enhance  the  transformer  ratio.  Several  techniques  are  currently  un¬ 
der  consideration  ranging  from  transverse-to-  longitudinal  phase  space  exchanger  ,  ad  initio 
laser  pulse  shaping  ,  and  wakefield-induced  shaper  combined  with  ballistic  compression  . 
These  different  techniques,  foreseen  to  be  available  at  ASTA,  range  in  their  achieved  pre¬ 
cision  of  the  control  and  in  their  complexity.  These  scheme  will  be  throughly  investigated 
at  ASTA. 
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5  Bench  test  of  a  dielectric  structure 


In  slab  DLW’s  the  electromagnetic  modes  can  be  categorized  into  two  types  longitudinal 
section  electric  (LSE)  and  longitudinal  section  magnetic  (LSM)  corresponding  respectively 
to  the  boundary  conditions  Ey  —  0  and  Hy  =  0  at  the  vacuum-dielectric  interface  (the 
axis  configuration  is  shown  in  Fig.  28.  The  lowest-order  LSM  mode  is  typically  used 
for  producing  the  axial  electric  field  Ez  necessary  for  acceleration.  Fig.  27  presents  the 
dispersion  curve  k{u)  associated  to  the  LSM  modes. 


Figure  26:  Experimental  setup  used  to  characterize  DLW  structures  (photograph)  and 
geometry  of  the  DLW  under  consideration  (inset  schematic).  The  legend  is  as  follows  “CL”, 
“OAP”  and  “LN”  respectively  stand  for  “cylindrical  lens”,  “off-axis  parabolic  mirror”,  and 
“LiNbOs  crystal”.  The  red  lines  indicate  the  path  of  the  800-nm  4-mJ  IR  pulse  and  the 
yellow  lines  correspond  to  the  THz  beam. 

In  preparation  for  an  experiment  at  the  Advanced  Superconducting  Test  Accelerator 
(ASTA)  at  Fermilab  we  are  developing  a  THz  bench  test  for  slab  DLWs.  In  this  setup  a 
THz  pulse  is  used  to  probe  the  electromagnetic  properties  of  the  slab  DLW  using  a  near¬ 
field  THz  mapping  technique  similar  to  what  was  developed  in  Ref.  .  The  method  can  be 
used  to  obtained  spatial  distributions  of  the  modes  along  with  their  dispersion  relations. 

5.1  Experimental  Setup 

The  experimental  setup  to  be  used  to  characterize  the  DLW’s  is  depicted  in  Fig.  26.  The 
basic  idea  is  to  eventually  perform  a  near-field  mapping  of  a  guided  short  THz  pulse  at  the 
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Figure  27:  Dispersion  relation  associated  to  the  LSM  modes  in  a  slab-symmetric  DLW 
with  parameter  a  —  100  fim  and  b  —  120  iim.  The  red  line  correspond  to  kz  =  2nfc. 

input  and  the  output  of  the  DLW. 

A  4-mJ,  800-nm  amplified  laser  pulse  produced  by  a  commercial  laser  system  is  used  to 
produce  a  THz  pulse.  The  laser  is  obtained  from  regenerative  amplification  of  a  broadband 
5 X  ~  100  nm  pulse.  Several  method  to  produced  THz  were  attempted:  “conventional” 
optical  rectification  by  impinging  a  ZnTe  crystal  with  100  cut  and  phase-matched  optical 
rectification  using  a  LiNbCU  (LN)  wedged  crystal  .  Because  of  its  significantly  higher  yield 
and  relatively  easy  implementation  we  choose  the  latter  method  as  shown  in  Fig.  26  and 
use  a  stoichiometric  MgO(0.6%):LiNBO3  crystal  with  a  5  x  5  mm2  entrance-face  size. 

The  pulse  impinges  a  1200-lines/mm  grating  to  provide  the  necessary  wavefront  tilt.  A 
cylindrical  lens  focusing  in  the  horizontal  plane  (CL1,  /  =  15  cm)  provides  an  approximate 
image  of  the  grating  onto  the  LN  crystal.  A  second  cylindrical  lens  (CL2)  is  used  to  focus 
the  beam  vertically.  The  typical  spot  on  the  LN  crystal  face  is  wx  x  wy  ~  3  x  0.5  mm2. 
Upon  exiting  the  LN  crystal,  the  THz  beam  is  collimated  with  an  off-axis  parabolic  mirror 
(OAP1)  and  focused  inside  the  DLW  structure  (OAP2).  The  transmitted  THz  beam  is 
then  collimated  (using  OAP3  and  4)  and  sent  to  a  Michelson  interferometer  equipped  with 
a  Helium-cooled  InSb  bolometer.  In  a  later  stage  we  plan  to  have  an  electro-optical  imaging 
diagnostics. 

5.2  Simulation  and  Analysis  of  a  THz  pulse  propagation  in  a 
SLAB  DLW 

To  assess  the  performance  of  the  described  experimental  setup  we  modeled  the  propagation 
of  a  THz  pulse  through  the  DLW  and  use  these  numerical  results  to  develop  analysis 
tools.  The  simulation  of  a  THz  pulse  through  the  DLW  were  performed  using  the  program 
VORPAL  .  Vorpal  uses  the  conformal  finite  difference  time  domain  (CFDTD)  method  to 
solve  Maxwell’s  equations.  The  geometry  of  the  problem  simulated  with  VORPAL  appears 
in  Fig.  28. 

For  the  results  presented  below,  a  THz  pulse  is  launched  from  the  (x  =  0,y,z)  plane 
using  a  time  dependent  current  density 

Jz(t)  =  Aotsm((Jit)  sin(cc;2 1)  sin^ t)e~^~to^2^r2\  (37) 

where  uoi  =  2i rfi  where  (/i,/2,/a)  =  (0.8, 1.0, 1.2)  THz,  to  =  2  ps,  and  the  full-width 
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Figure  28:  A  diagram  of  the  slab  from  the  x-transverse  direction  (a)  and  the  from  the 
z-longitdinal  direction  (b).  The  slab  is  composed  of  a  dielectric  coating  with  dielectric 
permitivity  er  surrounded  by  a  perfectly  conductive  boundary  (PCB).  In  the  VORPAL 
simulation,  we  inject  the  THz  pulse  using  a  current  J  to  drive  a  short  pulse  in  z.  Lastly,  to 
remove  reflections  and  to  artificially  produce  the  pulse  leaving  the  structure  we  implement 
a  perfectly  matched  layer  (PML).  We  record  the  signal  on  axis,  near  the  PML. 


half-max  (FWHM)  duration  is  r  =  1  ps.  Although  the  pulse  peaks  at  frequencies  /),  its 
relatively  large  bandwidth  Sf  ~  r_1  1  THz  allows  for  the  excitation  of  higher  frequencies 

modes.  The  electromagnetic  field  of  the  guided  pulse  are  recorded  on  a  two-dimensional 
grid-line  as  a  function  of  time  at  the  entrance  Eu(x,y  =  80  /im,z  =  0,£)  and  at  the 
exit  Eu(x,y  —  80  /im,  z  —  2  cm,  t)  (where  u  E  [x,y,z]  corresponds  to  the  three  field 
components). 


Figure  29:  Fourier  transform  of  the  pulse  at  the  entrance  of  the  structure  |  Ey{x  — 
0,y  =  80  /im,z  =  0,  f)\  (top)  and  a  Fourier  transform  of  the  pulse  at  the  end  of  the 
structure | Ey{x  =  0,  y  =  80  /im,  z  =  2  cm,  /) |  (bottom). 

The  Fourier  transform  of  the  pulse  before  and  after  propagating  through  the  structure 
provides  a  partial  understanding  of  the  propagation  of  the  modes  through  the  structure. 
The  surviving  frequencies  correspond  to  the  eigenmodes  of  the  structure  while  components 
which  have  faded  can  be  associated  to  the  evanescent  modes.  Figure  29  shows  the  Fourier 
transform  of  the  Ey  component  at  the  entrance  and  exit  of  the  DLW  structure.  The 
spatial  variations  of  the  frequency  content  within  the  pulses  is  indicative  of  the  mode’s 
spatial  distributions. 

We  import  the  signal  recorded  at  the  end  of  the  structure  into  another  program  for 
spatio-temporal  analysis.  For  each  grid  point  recorded  on  the  transverse  line  at  the  exit 
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of  the  structure,  we  use  a  moving  window  to  scan  through  the  signal  temporally.  At 
each  window  step,  we  perform  a  Fourier  transform  and  record  the  strength  of  the  Fourier 
component  of  a  selected  frequency.  To  reduce  any  unwanted  artifacts,  e.g.,  from  sharp 
edges  the  moving  windows  is  a  Hann  window 


w(n)  =  0.5 


1  —  cos 


27 in  \ 
N  -l) 


(38) 


for  an  individual  point  n,  of  N  sampled  points. 

In  Fig.  30  we  show  the  spatio-temporal  distribution  of  the  guided  pulse  associated  to 
the  filtered  frequencies  of  1.4  and  1.8  THz.  Using  this  method,  and  comparing  to  the 
dispersion  curve  above,  we  can  confirm  the  propagation  of  the  LSM2m  and  LSM3m  for 
the  1.4  THz  component  as  well  as  the  propagation  of  the  LSM2m,  LSM3m,  and  the  near 
cutoff  region  of  the  LSM4m  which  explains  its  tardiness  passing  through  the  structure.  In 
addition  to  providing  information  on  the  modes’  spatial  distribution,  the  technique  can 
also  be  used  to  deduce  the  group  velocity  of  the  modes. 


horrizontal  position  (mm)  horrizontal  position  (mm) 


Figure  30:  Spatio-temporal  distribution  of  the  guided  pulse  (shown  in  Fig.  29)  filtered  at 
1.4  THz  (top)  and  1.6  THz  (bottom). 

5.3  Status  and  Conclusion 

An  experimental  setup  to  support  the  characterization  of  DLW  structure  was  assembled 
and  first  results  on  the  THz  pulse  production  were  recently  measured.  In  this  paper  we 
presented  numerical  simulation  to  benchmark  our  analysis  technique.  Our  results  shows 
that  quantitive  information  on  the  DLW  mode  can  be  obtained.  The  experiment  is  being 
carried  used  the  laser  available  at  the  high  brightness  electron  source  laboratory  (HBESL) 
at  Fermilab  National  Accelerator  Laboratory. 

6  Education  &  Collaboration 

One  student,  Mr.  Francois  Lemery,  is  currently  being  supported  full  time  on  this  grant. 
Mr.  Lemery  joined  this  project  in  the  Fall  2010  and  has  been  working  toward  numerical 
modeling  of  dielectric  wakefield  accelerators.  Mr.  Lemery  has  been  an  active  researcher  and 
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presented  his  work  at  various  international  conferences.  Most  recently  he  gave  two  talks  on 
electron-bunch  shaping  techniques  and  optimal  current  shapes  for  enhanced  transformer 
ratio  at  the  Advanced  Accelerator  Concept  Workshop  2014.  Mr.  Lemery  secured  a  one- 
year  University  dissertation-completion  fellowship  to  finalize  the  work  carried  under  the 
DTRA  grant.  He  will  graduate  in  Spring  2015. 

In  addition,  another  student  not  sponsored  by  this  grant  (but  supervised  by  P.  Piot)  is 
partially  involved  on  this  project.  Mr.  Christopher  Prokop,  who  joined  the  beam  physics 
group  in  the  summer  2008,  obtained  his  Master  in  the  summer  2009.  Mr.  Prokop’s  Master 
pertained  to  high-fidelity  particle-in-cell  simulations  of  table-top  Terahertz  coherent  light 
sources  based  on  the  Smith-Purcell  effect.  Mr.  Prokop  extensively  worked  with  the  code 
VORPAL  from  Tech-X  and  brings  his  expertise  by  helping  us  setting  up  the  problems  to 
be  simulated  as  part  of  this  project. 

This  project  was  a  joint  venture  with  Tech-X  and  we  have  continuously  collaborated 
with  P.  Stoltz,  D.  Smithe,  B.  Cowan,  and  J.  Smith.  We  have  especially  been  meeting 
regularly  (in  average  monthly)  to  discuss  our  progress  and  problems.  This  collaboration 
has  been  extremely  fruitful  and  led  to  several  publications  and  contributing  to  our  students 
training. 

As  part  of  this  proposal  we  strengthen  our  relation  with  the  German’s  national  lab¬ 
oratory  Deutsches  Elektronen-Synchrotron  (DESY)  located  in  Hamburg,  Germany.  We 
had  initially  approached  potential  collaborator  to  perform  a  DWA  experiment  there.  We 
were  granted  36  hours  of  beam  time  on  the  Free-electron  LASer  in  Hamburg  (FLASH) 
facility  in  2011.  This  is  a  commensurable  contribution  to  our  project  (operating  costs  of 
the  facility  are  $  5-10k/hour).  We  conducted  an  experiment  to  demonstrate  we  could 
shape  the  current  profiles  of  electron  bunches  to  enhance  the  performance  of  wakefield 
accelerators.  This  phase-I  experiment  resulted  in  very  exciting  results  that  was  published 
in  Physical  Review  Letters.  Subsequently  P.  Piot  visited  DESY  to  discuss  a  possible  DWA 
experiment  at  FLASH.  The  FLASH  team  was  very  receptive  to  such  an  experiment  and 
offered  significant  support. 

Finally,  we  collaborated  with  Euclid  Lab.  LLC  and  worked  with  personnel  at  this  small- 
business  innovative  research  company  to  specify  and  obtain  dielectric-wakefield  structures. 

7  Publications,  and  presentation 

7.1  Talks 

1.  F.  Lemery,  "Passive  longitudinal-phase-space  tailoring  of  non-ultrarelativistic  beams 
with  dielectric-lined  waveguides",  16th  Advanced  Accelerator  Concepts  Workshop 
(AAC  2014),  San  Jose,  CA,  July  13-18,  2014  (working  group  5). 

2.  F.  Lemery,  "Alternative  Shapes  and  Shaping  Techniques  for  Enhanced  Transformer 
Ratios  in  Beam  Driven  Techniques",  16th  Advanced  Accelerator  Concepts  Workshop 
(AAC  2014),  San  Jose,  CA,  July  13-18,  2014  (working  group  2). 

3.  F.  Lemery,  "Flat-beam-driven  Dielectric-wakefield  acceleration  at  ASTA",  ASTA 
users’  meeting,  July  24,  2013,  Fermilab.  Batavia  IL. 
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4.  P.  Stoltz,  ’’Modeling  of  thermal  effects  in  dielectric  wakefield  accelerators”,  APS 
March  meeting,  Baltimore  MD  (2013) 

5.  D.  Mihalcea,  ”A  dielectric-wakelfield  Energy  doubler  for  the  Fermilab’s  ASTA  facil¬ 
ity”,  High-brightness  electron  beam,  Puerto  Rico  (HBEB13) 

6.  P.  Piot,  ’’Formation  of  Tailored  Electron  Bunches  in  Multi-Frequency  Linacs” ,  invited 
talk  WG5,  Advanced  Accelerator  Concept  2012. 

7.  P.  Piot,  "An  Energy  doubler  DLW-based  module  for  the  ASTA  facility  at  Fermilab”, 
contributed  talk  WG4,  Advanced  Accelerator  Concept  2012. 

8.  P.  Stoltz,  ’’Wakefield  Generation  in  Compact  Rectangular  Dielectric  Loaded  Struc¬ 
tures  Using  Flat  Beams,”  53rd  Annual  Meeting  of  the  APS  Division  of  Plasma 
Physics,  Nov  14-18,  2011,  Salt  Lake  City,  UT. 

9.  P.  Piot,  "Generation  of  Electron  Bunches  with  Tailored  Current  Profiles  at  FLASH”, 
Bechleuniger  Seminar,  08/01/2011,  DESY-Hamburg,  Germany. 

10.  P.  Piot,  ’’Electron  bunch  shaping  for  enhanced  transformer  ratio:  ideas  &  issues”, 
DWA-FEL  workshop,  05/11/2011,  Argonne  IL. 

11.  D.  Mihalcea,  "High- gradient  wakefield  acceleration  in  dielectric-loaded  structures", 
Accelerator  Seminar,  05/05/2011,  Fermilab,  Batavia,  IL 
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Current  Profiles  in  a  Dual-Frequency  Superconducting  Linear  Accelerator",  Phys. 
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Dielectric-Lined  Waveguides",  accepted  Phys.  Rev.  ST  Accel.  &  Beams  (2014). 

4.  P.  Stoltz,  et  ah,  "Coupled  thermo-electromagnetic  simulations  of  beam-driven  dielec¬ 
tric  accelerator",  in  preparation 

7.3  Conference  proceedings  papers: 

1.  F.  Lemery,  and  P.  Piot,  "Passive  longitudinal-phase-space  tailoring  of  non-ultrarelativistic 
beams  with  dielectric-lined  waveguides",  Proceedings  of  the  16th  Advanced  Acceler¬ 
ator  Concepts  Workshop  (AAC  2014),  San  Jose,  CA,  July  13-18,  2014  (in  press). 

2.  F.  Lemery  and  P.  Piot,  "Alternative  Shapes  and  Shaping  Techniques  for  Enhanced 
Transformer  Ratios  in  Beam  Driven  Techniques",  Proceedings  of  the  16th  Advanced 
Accelerator  Concepts  Workshop  (AAC  2014),  San  Jose,  CA,  July  13-18,  2014  (in 
press) . 


final  report  —  HDTRA-1-10-1-0051 


page  42  of  43 


3.  F.  Lemery,  D.  Mihalcea,  P.  Piot,  J.  Zhu,  "Planned  High-gradient  Flat-beam-driven 
Dielectric  Wakefield  Experiments  at  the  Fermilab’s  Advanced  Superconducting  Test 
Accelerator”,  Proceedings  of  the  2014  International  Accelerator  Conference  (IPAC14), 
Dresden,  Germany,  June  15  -  20,  2014  (in  press). 

4.  F.  Lemery,  and  P.  Piot,  "Temporal  Electron-bunch  Shaping  from  a  Photoinjector  for 
Advanced  Accelerator  Applications",  Proceedings  of  the  2014  International  Acceler¬ 
ator  Conference  (IPAC14),  Dresden,  Germany,  June  15  -  20,  2014  (in  press). 

5.  F.  Lemery  et  ah,  "Toward  a  dielectric- wakefield  energy  doubler  for  Fermilab’s  ad¬ 
vanced  superconducting  test  accelerator",  Proc.  IPAC13,  in  press  (2013). 

6.  F.  Lemery  et  ah,  "THz  bench  test  of  slab-symmetric  dielectric  waveguide  structures", 
Proc.  IPAC13,  in  press  (2013). 

7.  F.  Lemery  et  ah,  "Performances  of  VORPAL-GPU  for  slab  DLW  structure",  Proc. 
IPAC13,  in  press  (2013). 

8.  P.  Piot,  et  ah,  "Formation  of  Tailored  Electron  Bunches  in  Multi-Frequency  Linacs", 
Proc.  Advanced  Accelerator  Concept  workshop  (AAC12),  in  press  (2012). 

9.  F.  Lemery,  et  ah,  "Optimization  of  Drive-Bunch  Profile  for  Enhanced  Transformer 
Ration  in  Beam-Driven  Acceleration  Techniques",  Proc.  of  the  2012  International 
Particle  Accelerator  Conference  (IPAC12),  3012  (2012). 

10.  F.  Lemery,  et  ah,  "Progress  toward  a  high-transformer-ratio  dielectric- wakefield  ex¬ 
periment  at  FLASH",  Proc.  of  the  2012  International  Particle  Accelerator  Conference 
(IPAC12),  2167  (2012). 

11.  D.  Mihalcea  et  al.  "Wakefield  Generation  in  Compact  Rectangular  Dielectric-Loaded 
Structures  Using  Flat  Beams",  Proc.  of  the  2011  Particle  Accelerator  Conference 
(PAC11),  340  (2011). 

12.  F.  Lemery  et  al.  "Experimental  plans  to  explore  dielectric  wakefield  acceleration  in 
the  THz  regime  ",  Proc.  of  the  2011  International  Particle  Accelerator  Conference 
(IPAC11),  2781,  (2011). 

13.  P.  Piot,  et  al.,  "Generation  and  Characterization  of  Electron  Bunches  with  Ramped 
Current  Profiles 


final  report  —  HDTRA-1-10-1-0051 


page  43  of  43 


DISTRIBUTION  UST 
DTRA-TR-1 5-36 


DEPARTMENT  OF  DEFENSE 


DEFENSE  THREAT  REDUCTION 
AGENCY 

8725  JOHN  J.  KINGMAN  ROAD 
STOP  6201 

FORT  BELVOIR,  VA  22060 

ATTN:  D.  PETERSEN 

DEFENSE  TECHNICAL 
INFORMATION  CENTER 
8725  JOHN  J.  KINGMAN  ROAD, 
SUITE  0944 

FT.  BELVOIR,  VA  22060-6201 
ATTN:  DTIC/OCA 

DEPARTMENT  OF  DEFENSE 
CONTRACTORS 

QUANTERION  SOLUTIONS,  INC. 
1680  TEXAS  STREET,  SE 
KIRTLAND  AFB,  NM  87117-5669 
ATTN:  DTRIAC 


DL-1 


